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This  research  concerns  linear  ordinary  differential  equations  depending 
in  such  a  way  on  a  parameter  p  that  the  "limit"  differential  equation 
obtained  by  letting  p  tend  to  *  in  the  differential  equation  is  of  lower 

\  r 

order  than  the  original  one.  f  t  C 

v  • 

Adopting  a  term  customary  in  physics  we  used  the  name  boundary  layer 
problem  for  the  question:  What  happens  to  the  solution  of  a  boundary  value 
problem  of  such  a  differential  equation,  if  the  parameter  tends  to  00  in  this 
solution? 

We  gave  a  general  answer  to  this  question  for  the  differential  equation 
■^N(y)  +M(y)  »  0,  where  N(y)  and  M(y)  are  linear  differential 
expressions  of  order  n  and  m,  respectively  (n  >  m) ,  and  for  non- 
homogeneous  boundary  conditions  which  consist  in  prescribing  the  values  of 
derivatives  (but  not  of  linear  combinations  of  such  derivatives)  at  the  end¬ 
points.  The  question  whether  the  solution  of  such  a  boundary  value  problem 
converges  to  a  solution  of  the  limiting  differential  equations,  as  p  ♦ 
and  what  boundary  conditions  are  satisfied  by  the  limit  function  could  be 
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decided  by  an  easily  applicable  rule.  This  rule  showed,  among  other  things 
that  the  solution  converges  only,  if  the  prescribed  n  boundary  conditions 
are  not  too  unevenly  distributed  between  the  two  endpoints . 

If  the  order  ra  of  the  limiting  differential  equation  is  only  one  less 
than  the  order  n  of  the  original  differential  equation,  then  the  above 
mentioned  rule  could  be  extended  to  more  general  types  of  boundary  conditions 
and  also  to  non-homogeneous  differential  equations. 

Since  the  most  important  boundary  layer  problems  in  the  applications  are 
concerned  with  systems  of  differential  equations,  we  gave  a  simple  example  for 
the  mathematical  treatment  of  a  boundary  layer  problem  for  a  linear  system  of 
two  simultaneous  differential  equations. 

The  validity  of  the  general  rule  proved  in  this  research  was  seen  to  be 
restricted  by  the  assumption  that  the  coefficient  of  the  term  of  highest  order 
of  differentiation  in  M(y)  has  no  zeros  in  the  interval  of  integration.  In  a 
special  example  we  showed  that  interesting  results  can  be  obtained,  if  this 
assumption  is  dropped. 

The  theory  of  the  asymptotic  expansion  of  the  solutions  of  linear 
differential  equations  involving  a  parameter,  developed  by  G.  D.  Birkhoff, 
Noaillon,  Tamarkin,  Trjitzinsky  and  others  proved  an  important  and  powerful 
tool  in  these  investigations. 
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In  May  1980,  the  Mathematics  Research  Center  organized  a  successful 
advanced  Seminar  on  Singular  Perturbations  and  Asymptotics  in  honor  of  the 
retirement  of  a  colleague,  Wolfgang  R.  Wasow.  His  fundamental  research  is 
responsible  for  many  other  rapid  developments  in  this  field  since  1940,  and 
continues  to  play  a  vital  role  in  modern  theory  and  current  applications. 
Wasow' s  Ph.D.  dissertation  (N.Y.U.,  1941),  a  small  part  of  which  exists  in 
print  (On  the  asymptotic  solution  of  boundary  value  problems  for  ordinary 
differential  equations  containing  a  parameter,  J.  of  Mathematics  and  Physics 
12  (1944),  173-183),  represents  the  starting  point  of  this  Important  flourish 
of  modern  applicable  research. 

Following  suggestions  of  several  participants  MRC  is  printing  his  1941 
thesis  in  its  entirety  as  a  TSR  in  order  to  make  this  valuable  work  more 
widely  available.  Readers  will  note  that  the  name  "singular  perturbations" 
(which  was  only  coined  several  years  later  by  K.  0.  Friedrichs  or  w.  Wasow  or 
possibly  jointly,  but  neither  is  now  able  to  recall  the  details)  does  not 
appear  anywhere  explicitly! 
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ON  BOUNDARY  LAYER  PROBLEMS  IN  THE  THEORY  OF 
ORDINARY  DIFFERENTIAL  EQUATIONS 

Wolfgang  R.  Wasow 
Introduction 

Many  problems  In  applied  mathematics  lead  to  questions  of  the  following  type: 

Given  Is  a  differential  equation  Involving  a  parameter 
p.  This  parameter  occurs  In  such  a  way  that  the 
’limiting*  differential  equation,  i.e.  the  differential 
equation  obtained  by  letting  p  tend  to  Infinity  In 
the  differential  equation,  is  of  lower  order  than  the 
original  one.  What  happens  then  to  the  solution  of  a 
boundary  value  problem  of  the  original  differential 
equation,  if  p  tends  to  infinity  In  that  solution? 

It  Is  by  no  means  obvious  -  and  not  even  always  true,  as  we  shall  see  -  that  the 
solution  of  such  a  boundary  value  problem  tends  to  a  solution  of  the  limiting  differential 
equation,  as  p  tends  to  Infinity,  But  even  when  this  is  the  case  the  question  arises 
what  are  the  boundary  conditions  satisfied  by  the  limiting  function.  As  a  solution  of  a 
differential  equation  of  lower  order  than  the  original  one  it  cannot,  in  general,  be 
expected  to  satisfy  all  the  boundary  contions  prescribed  in  the  original  problem. 

In  those  cases  in  which  the  solution  of  the  original  problem  converges  -  as  p  tends 
to  infinity  -  to  a  solution  of  the  limiting  differential  equation  which  does  no  longer 
satisfy  all  the  originally  prescribed  boundary  conditions,  the  solution  f  the  original 
problem  shows  a  peculiar  behavior  for  very  large  values  of  the  parameter  p.  Some  of  the 
derivatives  of  the  solution  will  assume  very  large  values  in  a  narrow  region  near  the 
boundary.  As  p  tends  to  infinity,  these  derivatives  will  tend  to  infinity  at  a  certain 
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part  of  the  boundary.  In  the  most  important  applications  of  phenomena  of  this  type  the 
first  derivative  of  the  solution  -  and,  of  course,  all  the  higher  derivatives  -  diverge  at 
parts  of  the  boundary,  as  p  tends  to  infinity. 

In  the  physical  interpretations  this  means  the  occurrence  of  "Boundary  layers"  in 
which  the  quantity  to  be  investigated  increases  or  decreases  very  rapidly  with  the  distance 
from  the  boundary,  when  some  physical  constant  is  large.  We  shall  use  the  name  Boundary 
Layer  Problems,  in  a  more  general  sense,  for  all  related  mathematical  problems. 

The  most  famous  problem  of  this  type  is  the  relationship  between  the  theories  of 
viscous  and  ideal  liquids.  An  interesting  boundary  layer  problem  for  a  system  of  two  non¬ 
linear  ordinary  differential  equations  has  been  investigated  recently  by  X.  Friedrichs  and 
J.  J.  Stoker  in  a  paper  on  the  buckling  of  elastic  plates,  [9], 

The  majority  of  the  applications  lead  to  non-linear  partial  differential  equations 
which  are  so  complicated  that  a  complete  mathematical  treatment  has  not  yet  been  attempted. 
But  even  the  boundary  layer  problem  for  linear  ordinary  differential  equations,  a  problem 
interesting  from  the  mathematical  as  well  as  from  the  physical  point  of  view,  has  as  yet 
been  hardly  investigated.  The  only  paper  known  to  the  author  of  this  investigation,  on 
this  problem  for  ordinary  differential  equations  is  the  article  [7]  by  Erich  Rothe,  in 
which  the  problem  is  solved  for  a  very  special  linear  differential  equation  of  the  second 
order  with  constant  coefficients. 

In  Chapter  I  of  the  present  paper  we  discuss  the  boundary  layer  problem  for  linear 
homogeneous  differential  equations  depending  linearly  on  the  parameter,  and  for  non- 
homogeneous  boundary  conditions.  The  result  of  this  part  can  be  expressed  by  a  simple  and 
easily  applicable  rule  which  determines  immediately,  for  a  given  problem  of  this  type, 
whether  the  solution  converges  and  what  boundary  conditions  are  satisfied  by  the  limiting 
function. 

In  Chapter  II  we  investigate  more  thoroughly  the  case  where  the  order  of  the  limiting 
differential  equation  is  lower  by  one  than  the  order  of  the  original  differential  equation. 
In  this  case  the  statement  of  the  general  rule  of  the  Main  Theorem  in  Chapter  I  can  be 
formulated  so  as  to  include  more  general  boundary  conditions  than  those  assumed  in  Chapter 


I.  The  problem  is  then  solved  -  at  leact  for  a  drop  of  one  in  the  order  of  the 
differential  equation  -  for  the  non-hoe  jeneous  equation.  Finally,  the  usual  method  of 
treating  boundary  layer  problems,  consisting  of  a  transformation  of  the  independent 
variable  and  an  appropriate  modification  of  the  boundary  conditions,  is  shown  to  be 
justified  in  this  case.  The  method  is  sometimes  used  in  more  complicated  problems  without 
mathematical  justification. 

In  Chapter  III  some  special  examples  of  other  types  of  boundary  layer  problems  are 
discussed. 

The  methods  used  in  this  paper  are  based  on  the  theory  of  the  asymptotic  solution  of 
ordinary  differential  equations  involving  a  parameter.  This  theory  has  been  developed  by 
G.  D.  Birkhoff  [1],  Noaillon  [2],  Tamar kin  [3],  [4],  Trjitzinsky  [6],  and  others.  In  the 
Appendix  we  give  a  short  outline  of  the  results  of  this  theory  as  far  as  they  are  used  in 
this  investigation. 

I  am  deeply  indebted  to  the  Professors  R.  Courant  and  K.  0.  Friedrichs  'hose  help  and 
encouragement  played  a  major  part  in  the  preparation  of  this  thesis.  The  original 
suggestion  for  this  investigation  came  from  Prof.  Friedrichs,  and  his  active  interest  in 
the  progress  of  the  work  has  been  of  the  utmost  value. 
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Chapter  I 


THE  BOUNDARY  LAYER  PROBLEM  FOR  THE  DIFFERENTIAL 

EQUATION  -  N(y)  +  M(y)  -  0 
P 

$1.  Statement  of  the  Problem 
1.  He  consider  the  ordinary  linear  differential  equation 

L(y,p)  «  0  , 

where  the  linear  differential  expression  L(y,p)  is  of  the  form 

L(y,p>  =  -  N(y)  +  M(y> 

P 


with 


(101) 


(102) 


n  (n-v) 

N(y)  =  l  a  (x)  y  (x)  (103) 

vO  V 

m  (m-y) 

M(y)  =  l  b  (X)  y  (x)  .  (104) 

P-0  U 

x  is  a  real  variable  and  p  a  positive  parameter.  He  assume  that  the  coefficients 

a  (x)  and  b  (x)  admit  at  least  n  bounded  derivatives  in  the  interval 
\)  u 

a  «  *  <  B  . 

If  the  order  n  of  the  differential  expression  N(y)  is  greater  than  the  order  m 
of  the  differential  expression  M(y),  then  the  differential  equation  (101)  gives  rise  to  a 
boundary  layer  problem  for  the  "limiting"  differential  equation,  i.e.  the  differential 
equation  obtained  by  letting  p  tend  to  infinity  in  the  original  differential  equation 
L(y,p)  -  0.  For,  this  limiting  differential  equation  is 

M(y)  -  0  ,  (105) 

and  this  differential  equation  is  of  lower  order  than  (101),  if 

n  >  m  (106) 

He  shall  also  assume  that 

m  >  0  .  (106a) 
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Most  of  our  results  reaaln  valid  for  a  ■  0.  But  at  some  points  the  Inclusion  of  the 
case  m  -  0  would  make  the  statement  of  the  result  rather  Involved.  It  seemed  therefore 
preferable  to  exclude  this  case  from  the  Main  Theorem. 

Together  with  the  differential  equation  (101)  we  prescribe  n  boundary  conditions  for 


the  function 

y(x).  The  boundary  conditions 

considered  in 

this  chapter 

are  of  the  form 

Li<y>  “  4i  ' 

(i  -  1,2 . 

n) 

(107) 

with  constant 

and  with 

(  <Xi‘ 

1 

y  <6) 

for  i  “  1,2, 

r  •  «  »  /I 

4 

Vy)  *  {  ,  , 

(108) 

* 

l  y  fa) 

for  i  “  r+1. 

r  e  •  »  ,11  , 

where  x  -  a  and  x  »  8  are  the  left  and  right  endpoints,  respectively,  of  the  interval 
under  consideration.* 

We  assume  that  the  boundary  conditions  are  arranged  in  such  a  way  that 

X.  >  X,  >  •••  >  1 
1  2  r 


Vi  >  Va  >  —  >  Tn  • 

This  arrangement  is  the  opposite  of  the  customary  one,  but  it  is  more  practical  in  our 

case.  Ml  the  numbers  1^  and  are,  of  course,  assumed  to  be  less  than  n.  r  is  the 

number  of  boundary  conditions  prescribed  at  the  right  endpoint.  The  number  of  boundary 

conditions  at  the  left  end  point  is  then  n  -  r. 

One  or  both  of  the  numbers  X  and  t  may  be  zero,  which  means  that  the  value  of 

r  n  J 

the  function  itself  is  prescribed  at  one  or  both  endpoints.  But  our  theory  applies  also  to 
cases  in  which  only  derivatives  of  the  function  are  prescribed  at  the  endpoints.  The 
boundary  conditions  (108)  contain  as  a  special  case  the  initial  value  problem.  We  have 
only  to  set  r  “  0,  or  r  “  n. 


SI 


*  dK 

we  shall  use  throughout  this  paper  the  notation  y*  ** ( x)  for  — y  , 

dx* 


We  make  further  the  assumption 


Sq(x)  9<  0,  for  ail  x  in  o  <  x  <,  6  ,  (109) 

which  makes  it  possible  for  us  to  set 

a0(x)  -  1  .  (110) 

without  loss  of  generality. 

A  very  essential  condition  for  the  validity  of  the  theory  that  follows  is  that  we  must 
have  also 

bQ(x)  f  0,  for  all  x  in  a  <  x  <  6  •  (111) 

It  is  easy  to  see  that  a  theory  of  boundary  layer  problems  which  does  not  assume  (111)  mist 
be  expected  to  be  much  more  complicated.  For  bQ ( x)  is  the  coefficient  of  the  first  term 
ii  the  limiting  differential  equation  (105).  Hence,  if  bQ(x)  has  zeros  in  a  <  x  <  B, 
these  zeros  will,  in  general,  be  singular  points  of  the  solutions  of  the  limiting 
differential  equation. 

To  these  assumptions  we  will  have  to  add  two  more  conditions  of  a  rather  essential 
nature.  Since  these  conditions  can  be  more  easily  formulated  in  connection  with  our  Main 
Theorem  we  postpone  their  statement  for  a  few  pages. 

In  general,  the  differential  equation  (101)  will  have  a  unique  solution  U(x,p) 
satisfying  the  boundary  conditions  (107).  U(x,p)  depends  on  the  value  of  the  parameter  p. 
We  will  be  able  to  give  a  general  rule  which  allows  us  to  decide,  for  a  given  problem, 
whether 

u( x )  -  lim  U(x,p)  (112) 

p+<» 

exists,  and  which  are  the  boundary  conditions  satisfied  by  u(x).  We  shall  see  also  that 
u(x>,  when  it  exists,  is  a  solution  of  the  limiting  differential  equation  M(y)  “  0. 

The  behavior  of  U(x, p),  as  p  tends  to  infinity,  will  be  seen  to  depend,  in 
general,  on  three  things  only: 

(a)  On  the  number  n  -  m,  i.e.  the  difference  between  the  orders  of  the 
original  and  the  limiting  differential  equation. 
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(b)  On  r ,  i.e.  on  the  way  in  which  the  n  boundary  conditions  are 
divided  between  the  two  end  points* 

(c)  On  the  sign  of  the  coefficient  bQ(x)* 

2.  There  are  a  great  many  different  possible  cases  for  our  boundary  layer  problem*  In 
9ome  cases  U(x,  p)  converges,  as  p  ♦  in  some  cases  it  diverges,  and  there  are  boom 
special  occurrences  that  are  not  covered  by  the  Main  Theorem*  This  accounts  for  the  fact 
that  the  Main  Theorem,  although  very  simple  to  apply,  is  somewhat  lengthy  to  formulate*  He 
precede  its  general  formulation  by  a  few  examples,  in  order  to  give,  without  proof,  an  idea 
of  the  variety  of  possible  occurrences*  In  the  convergent  cases  the  boundary  conditions 
satisfied  by  the  limit  function  u( x)  are  obtained  by  canceling  n  -  m  of  the  given 
boundary  conditions,  usually  taken  among  those  involving  higher  orders  of  differentiation* 
Example  1 . 

L(x,p)  z  —  y"'  +  x3y  +  2  x  y"  »  0 
P 

with  the  boundary  conditions 

y"  (  a)  -  l3  y"(“)  “  ^ 

y’lfl)  *  t2  * 

If,  e.g.,  a  ■  1,  8  “  2,  then  bQ(x)  >0  in  a  <  x  <  8  and  the  solution  U(x,p)  of  the 
problem  tends  to  the  solution  of  the  differential  equation 

M(y)  s  x3y"  ♦  2  x  y  -  0 

satisfying  the  boundary  conditions 

y"  (  8)  -  t, 
y’<8)  -  *2  , 

which  are  obtained  by  canceling  the  boundary  condition  given  at  x  “  a*  If  a  ”  -2, 

8  =  -1,  then  bQ(x)  <  0  and  U(x,p)  tends  to  the  solution  of  M(y)  •  0  with  the 
boundary  conditions 

-7- 


irri'-  "*  TT  ~‘“1 "  •  *  ^ 


I 


1 

* 


y-(a)  -  ^ 


yMB)  -  i. 


obtained  by  canceling  the  first  boundary  conditions  at  x  «  0. 

But  if  a  “  -I,  B  “  1 «  the  condition  (111)  is  no  longer  satisfied  and  our  Main 
Theorem  does  not  apply • 


Example  2. 


L(x,  p)  =  (y*4*  +  cos  x»y*3*  )  xy"  +  xy  -  0 


with  the  boundary  conditions 


y"<a)  -  i. 


y'(o)  -  i. 


If  a  “  -2,  B  “  -1.  then  b- (x)  <  0  and  0( x,  p)  converges  to  the  solution  of  the 


differential  equation 


satisfying  the  boundary  conditions 


y*  +  y  ■  0 


y*  (  o)  “  tj 
y  (o)  - 


which  are  obtained  by  canceling  the  first  boundary  condition  at  each  endpoint.  If  a  >  0, 
B  >  0  then  bQ(x)  >  0  in  a  <  x  <  8,  and  U(x,  p)  tends  in  general,  to  the  solution  of 
y*  +  y  »  0  with  the  boundary  conditiras 


because  in  this  case  the  general  theory  requires  the  canceling  of  the  two  boundary 
conditions  involving  the  highest  order  of  differentiation. 


a “  2*  ,  B “  , 

then  we  have  an  exceptional  case.  Because  then  there  is  no  solution  of  y"  ♦  y  ■  0 
satisfying  the  boundary  conditions  y(  cO  -  i^,  y(  B)  “  unless  ^  »  0.  Again,  our 

Main  Theorem  does  not  cover  these  special  values  of  a  and  B. 


Example  3. 


My,  p)  =  1  y<4)  -  xy'  -  y  -  0 


•>  '  f  iwyar**--  ^  **#««** 


with  the  boundary  conditions 


y’  "(a)  -  Jt4  y'  "(8)  -  Jt, 

y'  (B)  -  tj 
y  ( 6)  ■  *3  • 

If  a  <  0,  0  <  0,  then  (x)  >  0  and  U(x,p)  tends  to  a  solution  of 

xy*  -  y  -  0 

with  ths  boundary  condition  y(  0)  -  ij ,  because  the  Main  Theorea  requires  the  canceling  of 
two  boundary  conditions  at  the  right  endpoint  and  of  one  boundary  condition  at  the  left 
endpoint*  But  if  a  >  0,  0  >  0,  i.e.  bQ(x)  <  0,  then  U(x,p)  does  not  converge  at  all* 


(101) 

(107) 

(102) 

(103) 

004) 

(108) 
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Ha  make  the  following  assumptions] 


1°.  x  is  a  real  variable. 

2°.  p  is  a  real  positive  parameter. 

3°.  The  real  functions  a  (x)  and  b  (x)  have  at  least 

v  u 

n  bounded  derivatives  in  the  interval 

a  <  x  <  8  « 


4°.  n  >  m  >  0  . 

5°.  Bq(x)  ■  1  . 

6°.  bg(x)  y  0,  for  all  x  in  a  <  x  <  8  * 

7°.  n  >  X1  >  Ij  >•••>  \r  >  0 

n  >  Tr+1  >  Tr+2  >,‘*>  Tn  >  0  * 

Then  the  behavior  of  U(x,p),  as  p  tends  to  infinity,  can  be 
found  by  the  following  procedure: 

First  Step.  Find  the  remainder  s  of  the  division  of  n  -  m  by  4. 
Second  Step.  Find,  in  the  table  on  the  next  page,  the  values  of  the 
numbers  q  and  p  for  the  differential  equation  under  consideration. 


©  If  a  -  1,  bQ  >  O' 
or 

a  ■  3#  bQ  <  0- 

(il)  If  3  »  1,  b0  <  o' 
or 

s  -  3,  bQ  >  0 

(ill)  If  s  »  0,  b0  >  0^ 
or 

s  =  2,  b-  <  0  y 


n-mt1  n-m-1 

then  q  »  — - — ,  p  -  — - — 


. .  n-m-1  n-m+1 

then  q  »  — - — ,  p  -  — - — 


. «  n-m  n-ra 

then  q  »  — ,  p  -  — 


(113) 


(IV)  If  s  -  0,  bQ  <  01 


8  -  2,  bQ  >  0 


n-m- 2  n-m-2 

then  q  »  — - — ,  p  -  — - — 


■fl 


t 

Mi 

* 


Third  Step. 

A)  If  the  differential  equation  under  consideration  belongs  to 
one  of  the  cases  0  -  (iTl)  *  of  the  table  above ,  try  to  cancel  p 
of  the  boundary  conditions  at  the  point  x  “  B  and  q  of  the 
boundary  conditions  at  the  point  x  •  a,  going  in  each  group  of 
boundary  conditions  from  those  containing  higher  derivatives  to  those 
with  lower  derivatives .  This  is  only  possible,  if  there  are  enough 
boundary  conditions  on  either  side  to  be  canceled. 

B)  If  the  differential  equation  under  consideration  is  of  the 
type  IV,  proceed  first  as  under  A).  Prom  the  remaining  boundary 
conditions  cancel  then  those  two  which  contain  the  highest  order  of 
differentiation  without  regard  to  the  endpoint  at  which  they  are 
given.  It  can  happen  that  the  boundary  conditions  to  be  canceled 

in  application  of  this  last  rule  are  not  uniquely  determined,  because 
one  would  have  to  decide  between  two  boundary  conditions  of  the  same 
order  of  differentiation.  He  shall  call  this  last  occurrence  the  "Case 
of  Indetermination". 

Convergent  Case.  If  it  is  possible  to  apply  the  rule  of  the  Third 
Step  in  a  uniquely  determined  fashion,  then 

u(x)  -  lim  0(x,  p) 

P*** 

exists  and  is,  in  a  <  x  <  8,  that  solution  of  the  differential 


equation 


M(y)  -0 


which  satisfies  the  boundary  conditions  not  canceled  in  the  Third  Step 
of  this  rule,  provided  the  following  two  conditions  are  satisfied! 


He  use  the  circles  around  these  numbers 
distinguish  this  division  into  four  cases 
introduced  presently. 


,  writing  0 ,  (il)  ,  Qy)  ,  (w)  ,  in  ord 
from  another  division  into  two  cases  X, 


order  to 
X,  XX  to  be 


**~*er?~ 


8°.  If  the  boundary  conditions  not  canceled  in  tha  Third  Btap 
of  this  rula  ara  raplacad  by  tha  corresponding  homogeneous 
boundary  conditions,  than  tha  problem  determined  by  thaaa 
boundary  conditions  and  tha  diffarantial  aquation  M<y>  •  0 
haa  only  tha  solution  y(  x)  =0. 

8°.  Ho  two  of  tha  boundary  conditiona  cancalad  in  tha  Third 
Stap  at  x  «  a  have  ordara  of  diffarantiation  that  ara 
congruent  modulo  n  -  a,  and  tha  saaia  is  true  for  tha 
boundary  conditions  at  tha  right  and  point. 

Divergent  Casas.  If  tha  rula  of  tha  Wiird  Stap  cannot  be  applied 
because  at  one  of  the  endpoints  there  ara  not  enough  boundary  conditions 
to  be  canceled,  than  U(x,p)  will,  in  general,  not  converge,  as  p  ♦  «. 
The  proof  for  tha  divergence  given  in  this  investigation  la  valid  only 
under  two  assumptions  analogous  to  8°  and  9°,  which  for  their  formula 
tion  require  an  additional  remark: 

Fourth  Step.  If  the  rule  of  the  Third  Step  cannot  be  applied  because 
the  number  of  the  boundary  conditions  at  one  endpoint  is  smaller  than 
the  boundary  conditions  that  would  have  to  be  canceled,  then  cancel  all 
the  boundary  conditiona  on  this  side  and  so  many  boundary  conditions  on 
the  other  side  (going,  as  before,  from  higher  to  lower  order  of 
differentiation)  that  m  uncanceled  boundary  conditiona  remain.  Then 
we  make  the  assumptions: 

8°' .  If  the  boundary  conditions  not  cancelled  in  the  fourth  step 
of  this  rule  are  replaced  by  the  corresponding  homogeneous  boundary  con¬ 
ditions,  then  the  problem  determined  by  these  boundary  conditions  and  the 
differential  equation  M(y)  -  0  has  only  the  solution  y(x>  =  0. 

9  .  Ho  two  of  the  boundary  conditions  canceled  in  the  fourth  step 

at  x  “  a  have  orders  of  dlffsrantlation  that  are  congruent  modulo  n  -  m, 
and  the  same  is  true  for  the  boundary  conditions  at  the  right  endpoint. 


12 


T 


'MKfil 


i 
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Conclusion  In  the  divergent  case. 


1 )  If  the  rule  of  the  third  step  cannot  be  applied  because  at  one  of 

the  endpoints  there  are  not  enough  boundary  conditions  to  be  canceled,  and 

If  conditions  8°  and  9°  are  satisfied,  then 

lim  U(x,p)  -  ±«,  for  all  x  in  a  <  x  <  6  . 

P-*- 

2)  If  the  rule  of  the  third  step  cannot  be  applied  because  of  indeter¬ 
mination,  and  If  8°  and  9°  are  satisfied  for  each  of  the  two  possible 

ways  of  applying  the  cancellation  rule,  then  U(x,p)  does  not  converge,  as 
p  ♦  -  but  r sixains  bounded. 


4,  The  reader  is  advised  to  check  the  examples  given  in  fl  in  the  light  of  the  Main 
Theorem,  in  example  2,  in  particular,  we  discussed  a  case  in  which  assumption  8°  was  not 
satisfied.  We  now  give  an  example  where  assumption  9°  is  not  satisfied: 


Example  4.  n  ■  5,  m  «  2,  bn  >  0. 


y*  "  ( a)  -  t. 


y*  (a)  -  lA 


yI4>(8>  -  t. 


y‘  <&>  -  t. 


y  <o>  -  *5 

Here  n-m  ■  3,  hence  s  -  3.  From  table  (113)  we  find  q  ■  1,  p  «  2. 

The  two  boundary  conditions  that  are  to  be  canceled  at  x  »  8  have  the  orders  of 
differentiation  4  and  1.  But  4=1  (mod  n-m),  in  this  case.  This  means,  assumption  9°  is 
not  satisfied,  and  the  Main  Theorem  does  not  apply.  However,  if  bQ  <  0,  then  9°  is 
satisfied,  and  we  can  be  sure  of  the  convergence  of  U(x,p). 

It  is  an  open  question  whether  U(x,p)  can  converge  even  if  9°  is  not  satisfied.  It 
seems  unlikely  to  the  author  that  the  Main  Theorem  remains  valid  in  those  cases. 

The  next  example  is  of  the  type  which  we  have  called  the  case  of  indetermination. 


f  *•  **!f-*F- 


Example  5.  n  -  3,  ra  -  1,  bQ  >  0, 

y*  ( a)  -  i2  y<e)  -  t, 

y  («)  *  *3 

Here,  n-m  »  2  and  therefore  8-2.  Table  (113)  shows  that  this  is  the  case  ,  and 
that  p  -  q  -  0.  The  rules  of  the  Main  Theorem  would  require  the  canceling  of  the  two 
boundary  conditions  Involving  the  highest  derivatives.  This  cannot  be  done  in  a  uniquely 
determined  way,  since  y( 0)  -  just  as  well  as  y{  a)  -  might  be  canceled  In  addition 

to  y'(a)  -  .  Hence,  this  is  the  case  of  indetermination,  and  U(x,p)  does  not 

converge . 

The  rest  of  this  chapter  is  devoted  to  the  proof  of  the  Main  Theorem. 

§3.  Asymptotic  Solution  of  the  Differential  Equation  L(x,p)  -  0 

5.  As  pointed  out  in  the  introduction  the  principal  tool  of  our  proof  of  the  Main  Theorem 
is  the  theory  of  asymptotic  solution  of  differential  equations  involving  a  parameter.  He 
begin  by  defining  what  we  shall  understand  by  asymptotic  equality  in  this  investigation. 
Definition:  The  functions  f(x,p)  and  g(x,p>  are  said  to  be  asymptotically  equal  in  an 
interval  a  <  x  <  0,  if 


..  E(x,p) 

f  ( x ,  p )  -  g(x,  p)  +  - . 

P 

Here  a  >  0,  (but  not  necessarily  an  integer),  and  E(x,p)  is  a  function  such  that  there 
is  a  positive  real  number  R  so  that  |E(x,p)|  is  uniformly  bounded  for  o  <  x  <  0,  and 
p  >  R. 

If  a  function  f(x,p)  is  asymptotically  equal  to  a  function  F(x)  independent  of 
p,  we  shall  write 

f(x,p)  -  (F(x)J 
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Note  that  the  symbol  (F(x)J  does  not  describe  the  function  f(x,P)  uniquely, 
correct  to  conclude  from 

f^x.P)  -  [F(x)J 

and 

f2(*,P>  -  lr(x» 


that 


fj(x,P)  -  f2(x,p)  . 


It  is  not 


6.  Using  Noaillon's  method  the  following  theorem  can  be  proved. 

Theorem  1;  If  the  assumptions  1°  -  6°  of  the  Main  Theorem  are  satisfied,  then  the 
differential  equation  (101)  admits  a  complete  set  of  n  linearly  independent  solutions  of 
the  form 


Ou(*,P) 


f 


tn(x)]. 


(V 


1,2,... ,n-m) 


(u(x)  ]  ,  (V-  n-m+1,  n-m+2 ..... n ) 

^  v-n+m 


(114) 


(115) 


Here  we  are  using  the  following  abbreviations 

i/n__ 

1 )  a  —  |  p  | 

2)  ♦,(x),  ♦  (x),...,4  <x) 

I  i  n-ra 

n-m  values  of 


<-bQ(x)) 


1/ 


n-m 


arranged  in  such  a  way  that 


Re(^)  >  Re<$2>  • 


3) 


a  (Ob  (E)  -  b  (t) 

■ra  ' 


■) 


n(x) 


(116) 


(117) 


(118) 
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4 

4  i 

*  I 

fc*l 

* 


4)  The  function*  u^(x),  ( u  -  1,2, ...,■)  are  any  n  aolutiona  of  the  differential 


equation 


M(y)  -  0 


forming  a  complete  linearly  independent  system  of  such  solutions. 

The  equations  (114)  and  (115)  may  be  formally  differentiated  at  least  n-1  times. 


n^i>Cx,p) 


for  i  -  0, 1 , . . . , n-1 . 


l*v(x)*n(x)]  ,  (v  -  1,2, .  • .  ,n-m) 


*uv-n+m*  '  *v"  n-m+1 »  n-sM-2,...,n) 


7.  Remarks. 


As  roots  of  one  and  the  same  real  function,  the  complex  functions  *Mx)  are  of  a 
particularly  simple  structure.  If  bQ(x)  <  0,  then  the  *Mx>  ere  obtained  by 
multiplying  the  (n-m)-th  roots  of  unity  by  the  factor 


A  similar  relation  holds  when  b0(x)  >  0.  More  precisely:  Set 


,  if  bQ  <  0 


1/n_  »1 

lbQ(x)  |  en-m  ,  if  bQ  >  0 


Then  the  «^(x)  are,  -  in  different  order,  -  equal  to 

.  ,  .2  ,  n-m-1 

k,  ke,  ke  ,  ...»  ,  kt 


Furthermore,  sine*  the  *y(x)  art  the  root*  of  e  reel  nuaber,  the  sequence  of  the  ^(x) 
conalste  of  pelre  of  conjugate  nuaber a.  Hence,  In  (117)  the  •>"  end  the  signs 

alternate. 


§4.  Outline  of  the  Asymptotic  Solution  of  the 
Boundary  Value  Problem 

8.  He  try  to  repreeent  the  eolutlon  of  the  boundary  value  problem  defined  by  (101)  and 
(108)  in  the  fora 


n 

O(x.p)  -  l  c  (p)  0  (x,p)  .  (125) 

vl  v  v 

If  auch  a  eolutlon  exists,  then  the  function  cy(p)  are  solutions  of  the  system  of  linear 
equations 

n 

l  cv(p).L1(Oy)  -  tt,  (1  -  1,2, ...,n)  .  (126) 

Using  (108)  and  substituting  for  the  U(x,p)  the  expressions  (114)  we  find,  for 
v  *  1,2, ... ,n-a. 


w  - 


Xi  °"u  Xi 

a  e  (8)n(8)J.  (i  -  1,2, ...,r) 

^  ^  ,  ( v  ■  1,2, ... ,n-m) 

o  1  lFvX(<*)n(«)l.  (1  *  r*1,  r+2, . . . ,n) 


(127) 


where 


»„  -  /f  *U>dC  .  ( v  -  1,2,...,n-m)  . 

V  d  V 


For  v  «  n-mM, 


,n  we  have,  because  of  (115), 


W  ■  tWn-,i 


(128) 


(128) 


From  the  inequalities  (117)  it  follows  that  similar  inequalities  hold  for  the 
quantities  w^  (v“  1, 2, . . .  ,n-m) ,  i.e. 

Re(w„)  >  Re(w  )  >•••>  Re(w  )  .  (130) 

i  i  n-m 

In  order  to  find  the  c  (p)  of  (125)  from  (126)  we  have  to  calculate  the  determinant 


L^U,) 

L^Oj) 

•  •  •  ♦ 

LjfO,) 

•  •  • 

l2(u2) 

•  •  • 

•  •  •  •  L2(Un> 

•  •  • 

w 

•  •  • 

w 

....  Ln(U„) 

and  the  determinants  A^tp),  (v  “  1,2, 


,n),  obtained  by  replacing  the  v-th  column  of 


A( p)  by  the  column 


The  coefficients  cy(p)  are  then  given  by 


A  (p) 

c  (p)  “ .  y  ■  ■ 

v  A(  p) 


{ v  «  1,2, .. . ,n) 


(1.32) 


and,  if  we  substitute  (132)  in  (125)  we  obtain  the  form 


n  A  (p> 

0(x'p)  ”  \  MW  VX,p) 


for  the  solution  U(x,p)  of  our  boundary  value  problem. 
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9.  Our  aim  is  now  to  calculate,  with  the  help  of  tha  asymptotic  expressions  (114),  (IIS) 
and  (127),  (129)  the  asymptotic  value  of  the  right  number  of  (133).  The  first  and  most 
important  part  of  that  calculation  consists  in  finding  the  asymptotic  value  of  the 
determinant  dtp).  The  asymptotic  calculation  of  the  d^(  p)  does  not  offer  new 
difficulties. 


$5.  The  Asymptotic  Value  of  Mo) 

10.  As  a  consequence  of  (127)  and  (129)  all  the  terms  of  the  expansion  of  dtp)  are 
obviously  of  the  form 


1*1  d 


S  dV 


with  real  S  and  real  or  complex  V  and  K. 
Definition:  Two  expressions 


K,  -  t»t1l  a 


S,  dV, 


S  dv 

K2  -  (k2l  a  e 


with 


I*  0,  k2  ?  0 

and  real  S,  and  s2  will  be  said  to  be  of  equal  order  of  magnitude  if 

Retv,)  -  Re(Vj) 

and 


s,  -  S2 


If 


Re(V,)  >  Re(Vj) 


Re(V, )  -  Re(V2),  but  3,  >  Sj 

then  K,  is  said  to  be  of  higher  order  of  magnitude  than  K2,  and  vice  versa. 

If  R,  is  of  higher  order  of  magnitude  than  Kj,  then  we  can  obviously  write 

s  dv 

K1  ♦  k2  ■  ®  • 
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The  sum 

-  V,, 

l  [kj  oS  •  0  34) 

(I) 

of  all  the  t mtmM  of  highest  c  ler  in  A( p)  is  the  asymptotic  expression  of  A( p)  for 

large  p,  unless  all  the  V  are  alike  and 

Hi 

Ik  -  0  . 

L  m 

(i) 

In  this  latter  case  (134)  reduced  to 

.  .  s  vi 
(0]  o  • 

and  an  asymptotic  calculation  of  &{  p)  would  have  to  take  into  account  terras  of  lower 
order  of  magnitude  in  the  expansion  of  M p) ,  as  well  as  the  later  terms  in  the  asymptotic 
solutions  of  our  differential  equation.  The  exclusion  of  this  exceptional  case  from  our 
theory  will  compel  us  to  Introduce  the  conditions  8*  and  9*  of  the  Main  Theorem. 


11.  If  the  values  (127)  and  (129),  for  the  (U  are  substituted  into  the  expression 
(131)  and  A(p),  it  is  seen  that  the  last  m  columns  of  the  determinant  form  the  matrix 
of  n  rows  and  ra  columns 


(L1(u1)]  [L^  (Uj ) ]  •  •  ♦  [I^lu^)] 

(L^U,)]  tt,2(Uj)]  •  •  •  [Ljtu^)] 


|  [Ljjtu^)  H.„(u2))  •  •  •  [LjjtUjj)) 


(135) 


Ml  the  minors  of  this  matrix  have  an  order  of  magnitude  not  greater  than  that  of  1 . 

The  elements  of  the  first  n-m  columns  of  A(  p)  are  given  by  (127).  In  order  to 
find  the  asymptotic  value  of  A(  p)  we  expand  A(  p)  in  terms  of  its  n-m  first  columns 
and  investigate  the  order  of  magnitude  of  the  minors  in  this  expansion. 
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12.  Lemma  1 :  Let  D(p)  be  that  minor  of  the  determinant  M p)  which  is  formed  by  the 


first  n-m  columns  of  A(p,  and  by  its 


th  th 

11  '  S 


.  i 


th 


rows,  with 


i,  <  i2  <  •  •  •  <  in.m 


If  then  the  first  h  of  the  numbers  i  are  less  than  or  equal  to  r,  then  the  order  of 


magnitude  of  D(p)  is  not  greater  than  that  of 


h 

0  I  WN 

T  \m1 
a  e 


(136) 


where 


h  n-m 

*  -  I  \  *  1  \  •  f  137) 

v*  1  u  t#*h+ 1  p 

Proof i  From  (127)  and  (131)  we  see  that  each  of  the  first  h  rows  of  D(p)  contains  the 
factor  n(B)i  while  each  of  the  remaining  rows  contains  the  factor  n(a).  Hence,  we  can 
factor  in  D( p)  the  expression 


n-m-h 

n  (  S)  n  ( a) 


Furthermore,  we  see,  that  the  first  row  of  D( p)  contains  the  factor  o  ,  the  second 

Xi  Xi 

2 

row  the  factor  a  ,  etc.,  and  finally  the  h-th  row  the  factor  o  .  Similarly,  we  can 
factor  in  the  remaining  rows  of  D(p)  the  expressions 

Ti  Ti  Ti 

h+1  h+2  n-m 

o  ,0  ,  •  •  •  ,  o  , 

respectively.  Altogether  we  can  factor  in  D(p)  the  expression 


n-m-h 

h,  ..  ,  .  T 

n  ( 6)  n  ( a)  a  , 
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where  T  has  the  value  of  (137).  Then  D(p>  can  be  written  in  the  form 


1 

MS 


X. 

e  (8)1  •  • 


ow  1 

•  •  e  n  (6)1 

n-m 


n-m-h 

D(  p)  "  n  (  8)  n  («>  o 


[*,  n”m(a)  1 


"”m  ( a)  1 
n-m 


Now  we  expand  the  remaining  determinant  with  respect  to  the  minors  of  its  first  h  rows. 
Each  term  of  this  expansion  contains  an  exponential  factor,  and  none  of  these  exponential 
factors  is  of  greater  order  of  magnitude  than  the  one  originating  from  the  minor  formed  by 
the  first  n-m  rows  and  columns  of  the  determinant  of  (138).  For  the  exponential  factor 


in  this  term  is 


h 

o  l  * 


and  because  of  (130)  no  sum  of  h  w  's  has  a  greater  real  part  than  £  w^.  Hence,  the 

v  v»1 

order  of  magnitude  of  D(p)  is  not  greater  than  that  of 

h 

0  I  Wv 
T  \S>1 
0  e  . 

It  may  be  less,  for  we  have  to  take  into  account  the  possibility  that 


1  •*“ '  *  <fr+ 


V 


h 

o  l  » 

T  u*  1 
D(p)  »  fO  1  a  e 

Remark:  if  Re( wh+1 )  <  Re ( wh ) ,  than  there  la  only  one  term  of  maximal  order  of  magnitude 
In  the  expansion  of  the  determinant  of  (138)  with  respect  to  Its  first  h  rows.  But  If  we 
have  Re(wh+1)  -  Re(wh)  (compare  section  7  and  the  definition  of  the  w^,  formula  (128)), 
then 

h  h-1 

*«(  I  -v)  -  *•(  I  «v+  Vi)  ' 

\^1  ^1 

and  we  have  therefore  a  second  term  of  maximal  order  in  the  expansion  of  the  determinant  of 
(140)  with  respect  to  its  first  h  rows.  But,  clearly,  these  two  terms  cannot  cancel, 
since 

*h  *  Vi  • 

13.  The  question  of  finding  among  all  the  minors  of  the  first  n-m  columns  of  A(p)  the 
one  of  highest  order  reduces  now  to  the  two  questions: 

(a)  which  selection  of  n-m  rows  in  A( p)  leads  to  an  expression  (136)  of  highest 
order,  and, 

(b)  when  does  the  minor  corresponding  to  this  selection  actually  have  the  order 
indicated  by  (136). 

In  answering  the  first  question  several  cases  are  to  be  distinguished.  For  this 

distinction  the  number  of  v^ix)' s  which  have  a  positive  real  part  plays  an  essential 

role.  Let  us  call  this  number  p.  From  the  definition  of  the  V>  *s  it  is  clear  that  p 

v 

depends  on  the  sign  of  bg(x)  and  on  the  remainder  of  the  division  of  n-m  by  4.  A 
simple  calculation,  which  we  omit  here  shows  that  p  has  the  values  indicated  in  the  table 
(113). 

We  distinguish  the  following  cases: 

I.  Ret^^)  ft  0,  for  all  v  "  1,2,...,n-m 

A)  r-p  <  0 

B)  0  <  r-p  <  m 

C)  r-p  >  ra 
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II.  Ral^)  -  0,  for  some  v  . 

A)  r-p  <  0 

B)  0  <  r-p  <  m+2 

C)  r-p  >  o+2  . 

Remarks :  Remembering  the  definition  of  the  we  see  Immediately  that  case  II  of  (139) 

occurs  only  when 


n-o  =  0  (mod  4)  and  <  0 


(140) 


n-m  =  2  (mod  4)  and  b.  >  0 
0 


Case  II  Is  thus  seen  to  be  equivalent  with  the  case  (iv)  of  table  (113).  In  this  case 


© 


the 


>1 ways  exactly  two  functions  *  with  vanishing  real  parts.  In  the 


arrangement  of  (117)  these  are 


^(x)  and  i(x)  . 
p+1  p+2 


Casa  I  of  table  (139)  corresponds  to  the  cases  ^i)  -  of  table  (113).  in  these 

cases  the  number  q  of  table  (113)  is  equal  to  rt-m-p.  It  follows,  therefore,  that  the 
condition  IC  of  (139)  can  be  re-wrltten  In  the  form 

(n-r)  -  q  <  0  . 

In  case  of  (113),  which  we  have  seen  to  be  equivalent  to  case  II  of  (139),  we  see  from 

(113)  that 


q  »  n-m-p-2  , 

and  therefore  lie  is  also  equivalent  to  (141).  Hence,  we  can  say,  that  in  case  I  as  well 
as  in  case  II  of  (139) 

A)  is  the  case  in  which  the  canceling  rule  of  the  Main  Theorem  cannot  be 
followed  because  there  are  not  enough  boundary  conditions  at  x  ■  6. 

B)  is  the  case  in  which  the  canceling  rule  can  be  followed. 

C)  is  the  case  in  which  the  canceling  rule  cannot  be  followed,  because  there 
are  not  enough  boundary  conditions  at  x  ■  a. 
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r* 


14.  We  treat  the  case  I  of  (139)  first.  According  as  to  whether  we  are  in  the  case  IA, 


IB,  or 

IC, 

we  consider  then  the 

minor  D  (p),  D.  (p), 

a  d 

D0 ( p)  formed 

by  the 

n-m  first 

columns 

of 

the  determinant  A(p) 

and  by  the  rows 

1,2,  ee 

(case  A) , ( 142a) 

1,2,  .. 

.p> 

e  e  e  < 

r+1. 

r+2,  . 

e  ,  r+n-*-p  e 

(case  B) , ( 142b) 

r+2,  . 

so,  n  • 

(case  C) ,( 142c) 

respectively.  Using  lemma  1  we  shall  show  that  in  each  ease  the  einor  thus  defined  is  -  in 
general  -  of  the  highest  possible  order  of  magnitude  among  all  the  minors  of  the  first  n-m 
columns  of  A(p). 

Case  A)  The  expression  (136)  has  the  highest  possible  order t  if  the  rows  of  the  minor 

h 

are  chosen  in  such  a  way  that  J  w  has  the  greatest  possible  real  part.  In  case  IA 

v*1  v 

this  means  that  we  have  to  choose  h  •  r.  For  h  can,  -  by  definition,  -  not  be  greater 

than  r  and,  on  the  other  hand,  all  w^  with  v  <  r  have  positive  real  parts  in 

consequence  of  the  condition  r-P<0*  In  order  to  make  T  in  (136)  as  great  as  possible  we 

n-m 

have  then  to  choose  for  the  remaining  n-m-r  rows  of  the  minor  those  for  which  J  T 

ir»h+i  u 

is  greatest.  Since  the  are  arranged  in  order  of  decreasing  size,  this  is  the  case,  if 

we  choose  the  rows  r+1,  r+2,  ....  n-m.  This  is  exactly  what  we  have  done  in  (142a). 

h 

Case  B)  Here,  taking  h  -  r  would  not  make  the  real  part  of  J  w  a  maximum 

v-1 

because  we  have  r  >  p  and  the  sume  would  therefore  include  wv'*  with  negative  real 
part.  Instead,  we  have  to  take  h  ■  p  Including  thus  all  the  w^'s  with  positive  real 
part  and  only  those.  (137)  shows  then  that  taking  the  first  p  rows  of  A( p)  gives  the 
greatest  contribution  to  the  exponent  T.  For  the  remaining  n-m-p  rows  we  take  the  rows 
r+1,  r+2,  ...,  r+n-m-p  in  order  to  make  the  second  sum  of  the  right  number  of  (137)  a 
maximum.  This  is  possible  because,  in  consequence  of  (139),  r+n-m-p  <  n,  in  this  case. 


"H1 


Case  C)  Here  we  reason  as  In  case  B.  But  since  r+n-a-p  >  n  in  this  case,  takinq 
the  rows  r+1,  r+2,  ...,  n  in  addition  to  the  rows  1,2,  ...,  p  would  not  be  enough  to 
have  n-a  rows  altogether.  We  oust  therefore  choose  h  »  r-a,  in  order  to  have  n-a 
rows,  and,  as  before,  we  see  that  T  is  greatest,  if  we  take  the  first  r-a  rows  of 
A(p). 


Coaparing  (142b)  with  the  third  step  of  the  rule  of  the  Main  Theorea  we  see  that  the 
numbers  of  (142b)  are  just  the  subscripts  of  the  boundary  conditions  canceled  in 
application  of  the  Main  Theorem  in  the  cases  (7)  -  ^IIlJ ,  when  the  canceling  is  possible. 
Similarly,  (142a)  and  (142c)  contains  the  subscripts  of  the  boundry  conditions  canceled  in 
the  fourth  step  of  the  Main  Theorem  in  the  cases  -  ^1^ .  The  reason  for  this  fact 
will  appear  in  the  course  of  our  proof. 


IS.  The  considerations  of  section  14  are  not  sufficient  to  prove  that  the  ainors  D^tp), 

0.  (p),  D  (p),  respectively,  really  do  have  higher  order  than  all  the  others  (compare 

D  C 

section  12).  To  investigate  this  question  let  us  calculate  D  (p),  D.  ( p) ,  D  ( p) 

ADC 

explicitly.  Since  we  are  most  interested  in  case  IB,  which  we  shall  see  to  be  the 

convergent  case,  we  discuss  O.(p)  first. 

b 

In  order  to  calculate  D.  (p)  we  write  D.  (p)  in  the  fora  (138).  This  means  that  we 

D  D 

Interpret  the  numbers 

11'i2 . ^n-m 

of  (138)  as  being  the  numbers  (142b),  in  the  same  order.  Hence,  the  number  h  of  (138)  is 
in  this  case  equal  to  p.  From  the  definition  of  p  and  the  w^  it  follows  that 

S'fVl1  <  Re{wp>  * 

The  term  originating  from  the  minor  of  the  first  p  columns  in  the  expansion  of  the 
determinant  (138)  in  terms  of  its  first  p  rows  is  therefore  the  only  term  of  maximal 
order  of  magnitude  in  this  expansion.  We  obtain  therefore 

P 

a  I  w 

n-m-p  T.  **  v 

Dfe(p)  ■  (^(8)8  (a)  •  •  BfaJ  o  e  (143) 
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with 


and 


! 


x  ♦ 


V 


r+n-m-p 

l 


■t 

u 


(144) 


*  t  (B)  • 


•  •  *  (  B) 

P 


*,**(  B)  • 


•  •  *  MB) 
P 


(145) 


r*"— 

n-a 


v"1** 


„  r+n-a-p, 

Vi  ta) 


(146) 


(134)  is  valid  also  for  p  “  n-m  and  for  p  “  0,  If  M  define 
Bjj  -  1,  for  p  “  n-a.  p(x)  is  an  exponential  function,  hence 
D.  ( p)  has  therefore  the  order  of  magnitude  of 

D 


Aj,  -  1.  for  p  -  0, 
n(B)  »•  0,  n(u)  *  0. 


and 


if  Ab  ^  0  and  0. 

Lemma  2:  ^  0  if  and  only  if 

\i  i  A^,  (mod  n-a)>  (i,J  »  0,1,.. .,p)  (147) 

and  Bjj  f  0  if  and  only  if 

t  i  t.,  (mod  n-m)>  (lc,i  •  r+1,  . . . , r+n-m-p)  .  (148) 

n  1 


Proof:  The  statement  is  trivial,  as  far  as  is  concerned,  if  p  •  0.  For  p  >  0  we 

remember  from  $3,  section  7,  that,  with  the  notations  used  there,  the  *vfx)  are,  in 
different  order,  equal  to 


I 


k(x),  k( x)  e,  k( x)  t  ,  •  ...  ,  k( x)  e1*"""'  .  (124) 

In  tha  ordar  of  (124)  the  •f ^  ara  rapraaantad  in  tha  complex  plana  by  a  sequence  of 
successive  points  on  tha  circla  of  radius  |k|.  It  is  aasy  to  sea  that  tha  numbers  1^(0), 
¥>2(8)»  .••  i  v>p<  6)  ara  than,  in  different  order,  equal  to  tha  numbers 

k(0)et+1,  k<0)et4,2,  •  •  •  ,  ktB)®**-*  (149) 

where  t  is  a  certain  integer  which  is  only  determined  modulo  n-m.  Substituting  the 
expressions  (149)  into  (147)  we  find 


\^i 

i  k(0) 


(tMU. 


(ttl)X 


(t+p)  X, 
e 


(ttp)  X 
c 


(ISO) 


Now  we  set 


Cv'  C  ' 

This  allows  us  to  write  (150)  in  the  form 


(v*  1,2, «  ♦»  ,p) 


\p1 


^  -  *  k(6) 


t+1 


t+p 

•  •  •  C- 


t+1  tfp 

5i  '  *  •  s 


(151) 


where  is  the  Vandermonde  determinant  of  ^ Sinca  tha 

Vandermonde  determinant  vanishes  if  and  only  if  two  of  its  rows  are  equal,  ia  zero  if 

and  only  if  some  of  the  i;  are  alike.  But  from  (151)  we  see  that 


means 


C 


1 


or,  because  of  the  definition  of  e  in  formula  (123), 

Xi  =  Xj  ,  (mod.  n-m) 

This  proves  the  part  of  lemma  2b  that  is  concerned  with  The  proof  for  is  exactly 

analogous,  and  is  therefore  left  to  the  reader. 

The  reader  will  readily  see  that  similar  results  hold  for  D  ( p)  and  D  ( p) .  The 

4  c 

only  real  difference  in  the  reasoning  comes  from  the  fact  that  in  these  cases  there  may  be 
two  terms  of  maximal  order  in  the  expansion  of  (140).  But  since  these  terms  cannot  cancel, 
this  does  not  essentially  affect  our  argument,  we  restrict  ourselves  to  stating  the 
results  in  these  cases i 


T  0  l  Wv 

D#(p)  “  ffta)  aa.  ^ 

r-m 

o  £  w 
T  ^  v 

D (p)  -  [Q  J  a  c  e 
c  c 

where 


r 

I 

X  + 

n-m 

l  Tu 

v*1 

V 

l^r+1  u 

r-m 

l 

X  + 

n 

I  T 

V*1 

V 

irr+1  “ 

(  153) 


(154) 


(155) 


(156) 


and  Qa  and  Qc  are  two  constants,  with  respect  to  which  the  following  two  lemmas  hold: 
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Lemma  2a  i  Oa  ^  0  If 


and 


£  X^,  mod  n-m,  (i,j  «  1,2, ...,r) 


(157) 


£  Tj,  nod  n-m,  (k,I  -  r+1. 


,n-m) 


Lemma  2c:  Q_  *<  0  if 
—  c 


\  £  \y 

mod  n-m. 

<i,j  -  1,2,.. 

. . ,  r-m) 

Tk  *  V 

mod  n-m. 

(k, i  -  r+1, . , 

>•  «n) 

(158) 


Remark:  From  these  lemmas  and  the  rmurk  at  the  end  of  section  14  we  recognize  that  the 
conditions  (147)  and  (148)  are  equivalent  with  the  assumptions  9*  of  the  Main  Theorem  while 
(157)  and  (1S8)  are  equivalent  with  9*'.  Note,  however,  that  our  reasoning  so  far  does  not 
cover  the  case  II  of  (139),  which  we  have  seen  to  be  the  same  as  the  case  fiv)  of  (113). 


16.  The  cofactor  of  the  minor  D^tp)  in  4( p)  is  that  minor  of  the  matrix  (135)  which  is 
formed  by  the  m  rows  of  A(  p)  not  contained  in  D^(p).  let 


*b 


Vi'V . 

•  *  *  •  Vi'V 

W . 

•  •  •  •  L  (u  ) 
r  m 

Lr+n-m-p+1*u1 * 

•  •  •  •  X* 

r+n-m-p+ 1 

VV . 

•  •  •  •  L  (u  ) 
n  w 

Dtp)  is  of  the  form 

D 

(u  > 


(159) 


±  iy  . 

We  now  introduce  the  assumption 

^  +  0,  in  case  IB  .  (160) 

This  condition  will  be  seen  later  to  be  equivalent  with  the  assumption  8*  of  the  Main 


Theorem  in  the  cases 


Then  the  term  t  Db(p>  ( y  is  the  term  of  greatest 


order  of  magnitude  in  the  expansion  of  &(  p)  in  terms  of  its  first  n-a  columns.  Hence. 
(143)  gives 

T.  0  \  *v 

A(p)  “  [Qb  •  ^1  a  .  /  in  case  IB  /  (161) 


where 


n-m-p 

2b  -  ±rf(e>  •  n  (a)  •  •  Bjj  ?  0 


(162) 


Similarly,  if  we  define  and  6c  as  the  limits  -  as  p  of  the  matrix  (135) 

formed  with  the  rows  of  D(p)  not  occurring  in  D  (p)  and  0  (p),  respectively,  and  if 

B  C 

we  introduce  the  assumptions 


6  y  0,  in  case  IA 
a 

6  y  0,  in  case  IB 
c 


then  we  have,  in  analogy  with  (161) 


T  0  l  Wv 

A(  p)  -  •  6#]  •  o  e  , 

r-m 

_  o  y  w 

T  V 

A(  p)  •  [Q  •6]»oC«V"  * 

c  c 


in  case  IA  , 


in  case  IC  . 


(163) 

(164) 


(165) 


(166) 


This  finishes  the  asymptotic  calculation  of  Alp)  in  the  case  I  of  (139). 


§6.  The  Asymptotic  Value  of  the  Solution  of  the  Boundary 
Value  Problem  in  the  Case  IB. 

17.  The  method  used  in  $5  for  the  calculation  of  Alp)  can  also  be  applied  to  the 
determinants  A^(  p).  (For  the  definition  of  Ay(p),  see  §4,  section  8). 

For  v  >  n-m,  A^(  p)  is  distinguished  from  Alp)  only  in  one  of  the  last  m  columns. 
In  this  case  all  the  considerations  of  |5  remain  valid  for  Ay( p) ■  if  the  determinants 
6  ,  6  are  replaced  by  the  determinants  5^,  4^,  4,^  obtained  by  replacing  the 
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t 


i-(n-m)th  colvatn  of  4#,  5^,  by  the  cor  responding 


!•••>  by 


*n-a*1  *pM  1r-«+1 

•  •  • 

*  •  • 

•  ■  • 

•  l  or  • 

r 

•  jl  • 

rtn-n-p+l 


and  analogous  formulas  hold  in  the  two  other  cases. 
Hence 


(v  >  n-m)i  case  IB)  ,  (167) 


and  similarly  in  the  two  other  cases. 


yp>  ry 

&  (P)  ■  ly 


18.  For  v  <  nr-m  the  column  of  A(  p)  that  must  be  replaced  by  the  i^'s  in  order  to 
obtain  A  Ji  p)  changes  the  structure  of  the  determinant  somewhat .  But  if  we  place  this 
vth  column  behind  all  the  others  (and  change  the  sign  of  the  determinant,  if  necessary) , 
then  we  obtain  a  determinant  very  similar  to  A( p) •  The  only  essential  difference  is 
that  a  has  to  be  replaced  by  a'  ”  ih-1  (and  therefore  n-a  by  n-m-1)  and  that  w^ 
must  be  omitted  from  the  sequence 

*1<*2'  •  •  *  'wn-a* 
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T 


Again  we  treat  the  case  IB  (the  convergent  case,  as  we  shall  see)  first,  and  since  we 

V  p) 

are  more  interested  in  -  —  than  in  A Ji  p)  itself,  we  state  the  result  in  the  following 
font  i 

Lemma  3:  In  case  IB,  l.e.,  if  0  <  r-p  <  m,  we  have 


-X  -ow 

o  P  e  Ut  >  v  <  p 


r+n-m-p  ,  . 

a  [  it  ] »  p  <  v  <  n-m  , 


where  the  constants  depend  on  the  values  of  the  J^'s  and  are  not  necessarily 

different  froa  zero. 


SI 


1.)  v  <  p.  In  this  case  the  w^  to  be  omitted  has  a  positive  real  part.  Hence,  the 
formulas  for  A(  p)  can  be  applied  to  A  J  p)  if  we  replace 

p  by  p’  -  p-1 
m  by  m'  •  m+1 

and  the  sequence 


by  the  sequence 

w'l'w2 . ^m-1 

identical  with  the  first  one  except  for  the  omission  of  the  term  since  we  have 

0  <  r-p1  <  m1  , 

the  formula  to  be  used  in  (161).  He  obtain 

P’ 

T  0  £  Wu> 

Vp)  "  {Qbv  *  W  *  0  bV  6  UPl 


I 


i 


p'  r+n-a'-p' 

T m  l  X  *  l  T 

<i^1  ipr+1 


Note  that  the  reasoning  used  for  the  proof  of  ^  0  does  no  longer  hold  for 

Is  a  determinant  depending  on  the  .  It  Is  defined  as  the  cofactor  of  the  leading  minor 

In  the  expansion  of  fi  ( p)  with  respect  to  its  first  n-m-1  columns-  4  is  a 

V  D  V 

determinant  of  m+1  rows- 

From  the  definition  of  p'  ,  a',  w'  it  follows  thsn  that 


r-. ! 


)  w*  •  >  w  -  w 

u  •  w  v 
upI  up  1 


p-1  rfn-»*p 

t  -  y  x  ♦  j  t 

crl  irr+1  w 


Comparing  these  formulas  with  the  expression  for  M p)  we  see  that 


VI  f°bv  ‘  *bv)  _X, 

4  “>> "  l«b  •  vJ° 


p .  v 


2.)  p  <  v  <  n-a.  Hare  we  reason  as  under  1-),  the  only  difference  being  that 


p'  •  p  and  therefore 


r  -•  - ! 

u*  1  ur  i 


p  r+n-m-p- 1 

T.  -  ;  X  ♦  l  i 

bu  L.  u  u 

df>  1  ^  r+1 


Hence,  In  this  case 


•>  4  j 


3 


It  is  easy  to  see  that  the  first  two  sums  of  the  right  member  of  (169)  tend  to  sero,  in 
a  <  x  <  0,  as  p,  and  therefore  a,  tend  to  infinity.  Tor, 


R«(  /  -  w  )  <0,  in  a  <  x  <  0,  for  v  < 


fa  <  0,  in  a  <  x  <  B>  for  v  >  p  , 


in  consequence  of  (117),  (128)  and  the  definition  of  p. 

Therefore 

u(x)  «  lim  U(x,p)  -  £  ^ ^  u  '  (170) 

p+”  1^1  W 

But  if  we  remember  the  definition  of  «b,  ^  and  «Mx),  as  given  in  (159),  in 

section  17,  and  in  theorem  1  (section  6),  respectively,  then  we  see  that  the  right  member 
of  (170)  is  just  the  solution  of  M(y)  -  0  satisfying  the  boundary  conditions  not  canceled 
in  application  of  the  Hain  Theorem,  and  condition  (160)  is  seen  to  be  equivalent  with 
assumption  8°  of  the  Main  Theorem  in  the  convergent  cases  (T)  -  ^II)  (table  (113)), 


20.  Remark :  Formula  (169)  is,  in  fact,  a  complete  asymptotic  solution  of  the  boundary 
value  problem.  It  might  be  used  for  a  more  detailed  description  of  the  boundary  layer 
phenomenon.  It  can  be  written  in  the  following  simpler  and  more  symmetric  form: 

_T  P  of  mac 

U(x,p)  -  u(x)  +  a  T  £  e  (it  n(x)J 

\^i  v 


i  n-ra  of  *  (C)d£ 
♦  o‘X  l  e  v  ( 

\pp+ 1 


(*vn(x)] 


I 


&  • 


“1 

is 


where  T  «  t  and  X  »  X  are  the  lowest  orders  of  differentiation  occurring  in  the 

r+n-m-p  p 

canceled  boundary  conditions  on  each  side,  An  easy  consequence  of  (169)  is,  e.g.,  the 
following  interesting 

Corollary  1  Under  the  assuept ions  of  theorem  2  the  derivatives 

(t  -1) 

U(a,P ) ,  U*  (<*,P),  <  •  .,U(a,p)  r+n_“"P 

converges  to  the  value  of  the  corresponding  derivatives  of  u(x)  at  x  “  a.  The  next 
derivative, 

n  A  __  ) 

r+n-m-p 


0  (a,p) 


(t  > 

r+n-m-p 


converges,  as  P  +  •,  but  in  general  the  limit  is  not  equal  to  u  (a)  .  All  the 

higher  derivatives  of  U(x,P)  tend  to  infinity  at  x  ■  a.  An  analogous  statement  can  be 

(XJ 

made  at  x  -  0,  with  u(B,P)  as  the  last  convergent  derivative. 

In  less  precise  language  we  may  express  the  statement  of  this  corollary  by  saying  that 
the  last  canceled  boundary  condition  at  each  end  point  determines  the  derivative  of 
U(x,P)  in  which  the  boundary  layer  occurs  at  that  endpoint. 

§7.  The  Proof  of  the  Divergence  in  the  Cases  IA  and  IC. 

2t.  We  are  now  going  to  show  that  in  the  two  remaining  cases  U(x,P)  tends  to  infinity  as 

P  ♦  To  that  end  it  is  sufficient  to  prove  that  one  of  the  terms  in  the  right  member  of 

(133)  tends  to  infinity.  For  each  of  the  first  n-m  terms  of  (133)  is  of  the  form 

Y  a  Wv(x> 

Kva  e 


Y  0  Wjtx) 

If  one  of  these  terms,  say  K^o  e  tends  to  Infinity,  the  whole  sum  can  remain 

hounded  only  if  the  sum  of  all  terms  of  the  same  order  of  magnitude  as  this  particular  term 
vanishes  identically.  This  would  require  that  at  least  one  other  term  has  the  same 
exponential  factor  as  Kj.  We  would  have,  e.g.. 
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I 


a 


! 


w^(x)  »  Wj(x)  . 


Now  all  tha  H  (x)  are  of  the  fora 


V*>  ■  Vv*  /a  V£)dC 


where  the  constant  V  originates  from  the  factor  .  while  (x  <f  (E)d£  is  the 

v  A  (p)  Ja  v 

contribution  of  Uy(x,p).  But  an  equation  like 


V1  +  /a  VC)d5  “  v2  +  *2<5)d5 


is  impossible,  since  no  two  <f  (x)  are  equal. 


22.  Proof  of  the  divergence  in  the  case  IA: 

Let  us  calculate  the  determinant  A  ( p) .  In  this  case  the  formulas  for  A( p)  can 

r+ 1 

be  applied,  if  we  replace 

p  by  p'  ■  p-1  (since  r+1  <  p)  , 
m  by  m'  -  m+1  , 

and  the  sequence 


w1'w2 . wn-m 


by  the  sequence 


w1'w2' . wn-m-t 


obtained  by  omitting  the  term  wf+  ^  from  the  sequence  of  the  w^.  As  we  have 


r-p'  <  m1 


the  formula  to  be  used  is  (165).  we  obtain 


j  u 

.  .  .  .  ,  a, r+1  ufI 

Vi10’  ■  tQa,r+1  *  4a,r+15° 


It  follows  that 


r  n-m' 

T  .  ,  “  [  X  +  l  T 

4ir+1  «o  w  11 


r 

I 

Uf*1 


r 

l 

w*i 


and 


a, r+1 


r  n-m-1 

l  X  *  It 

*1  “  r*n  " 


For 


W» 

A(  p) 


we  obtain 


r+1 


A  ( 


<p)  fta.r+l  6a,r+l1  ”V« 

<»  "I  fi.  5.  J° 


and  the  r+1  at  term  of  (133)  becomes  therefore 


w> 

A  (p) 


“a,  r+1 


a. r+1 


nC  it) 


-  T 


r+1 


(0«C 


Since  r+1  <  p,  the  exponential  factor  of  the  right  member  tends  to  Infinity  as  p  ♦  «*• 

If  we  can  prove  that  the  expression  in  brackets  does  not  vanish,  then  the  divergence  of 
U(x, p)  is  assured,  n(x)  does  not  vanish  in  a  <  x  <  6  (see  (11S))«  For  fi4)r+1  we 

can  prove  that  the  following  lsmaa  is  true. 

Lemma  it  Q#  r+1  ?  0,  if  assumption  9*'  of  tho  Main  Theorem  is  satisfied. 

Proof:  An  almost  literal  repetition  of  the  arguments  of  $5,  applied  to  A^^ ( p)  instead 
of  to  A(p),  which  we  shall  omit,  shows  that  is  different  from  sero,  if  the 

following  determinants  do  not  vanish: 


39' 


X. 

*,(8) 


and 


a, r+1 


* 


(6) 


a, r+1 


Vi  Vi 

*  A  a) . *  (a) 

r+2  n-m 

Tn-*-1  Vm-1 

*  ,<a) . *  (a) 

r+2  n-m 


To  these  determinants  th«  reasoning  of  luma  2b  can  be  applied,  leading 

6  . ,  Is  a  determinant  one  column  of  which  is  formed  of  m+1  of 

•#rri 

1,.  8  . .  will  be  cero  for  certain  exceptional  vlues  of  these  t. . 

assumption  that  the  have  these  exceptional  values  would  not  be  sufficient  to  guarantee 
the  convergence  of  U(x,p),  for  the  r+1  at  term  of  (133)  will,  in  general,  not  be  the 
only  one  that  has  an  exponential  factor  tending  to  infinity*  We  shall  omit  the  not 
difficult  problem  of  finding  sufficient  conditions  for  the  convergence.  We  will  assume 
instead,  without  mentioning  it  each  time,  that  the  do  not  have  these  exceptional 

values. 

From  the  preceding  considerations  it  follows  that  U(x,p)  is  in  fact  divergent  in 
case  IA,  provided  that  the  assumptions  of  ths  Main  Theorem  are  satisfied. 


exactly  to  lemma  4. 
the  n  numbers 
But  even  the 


23.  The  proof  for  the  divergence  of  U(x,p)  in  the  case  ZC  is  analogous  to  the  proof  for 
the  case  ZA,  if  the  (r-m)-th  term  of  the  sum  in  (133)  is  considered,  instead  of  the 
(r+1 1st,  on  which  the  proof  in  the  case  ZA  was  based.  One  obtains  ths  expression 
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Ar-,(P> 

TIFT 


-P 


-S^ 


*  n(j.  )J 


-A  _  0{P  * 


(E)d£  -  w  } 
ot  r-m  r-m 


c  c 

from  which  the  divergence  of  U(x,P)  follows  as  in  the  case  IX,  since  Re(^r_m)  <  0. 


$8.  The  Case  II. 

24.  The  case  II,  i.e.  the  case  when  Reti^p^)  “  Re^p+2^  ”  °»  requires  a  special 

discussion,  because  in  our  reasoning  in  case  I  we  assumed  repeatedly  -  notably  in  section 

14  -  that  every  1*^(0^),  with  i  <  r,  is  either  of  greater  or  of  lower  order  of  magnitude 

than  all  L^lUy)  with  i  >  r.  But  this  is  no  longer  true  in  case  II,  for  the  Lt<Dp+i) 

°w  i  ow 

and  I‘1(Op+2>»  because  the  exponentials  e  and  e  ™  have  always  the  absolute 

value  1.  We  must  therefore  modify  our  considerations  for  case  II,  from  section  13  onward. 

We  remind  the  reader  of  the  remark  made  in  section  13,  to  the  effect  that  case  II  of 
(139)  is  equivalent  with  case  (iv)  of  (113). 


25.  The  case  IIX. 

Re-reading  sections  14-22  one  sees  that  no  modification  of  the  proofs  for  case  IX  is 

necessary  to  obtain  the  proofs  for  case  IIX.  Xll  arguments  remain  literally  the  same.  The 

Ow  i  «w 

reason  is  that,  for  r  <  p,  e  **  and  e  ***  do  not  appear  at  all  in  the  asymptotic 
expressions  for  A(p)  or  (P) . 


26.  The  case  IIB. 

This  is  no  longer  true  for  case  IIB,  defined  by  the  inequalities  0  5  r-p  <  m+2.  In 

order  to  find  the  minor  of  maximal  order  of  magnitude  among  the  minors  of  the  first  n-m 

columns  of  A(p)  we  go  back  to  lemma  1,  section  12.  Whether  the  number  h  in  (136)  is 

h 

chosen  equal  to  p,  p+1,  p+2,  does  not  have  any  influence  on  the  real  part  of  I  w  , 

h  v-1 

but  any  other  value  of  h  leads  to  a  \  w  with  a  smaller  real  part.  Which  of  these 

v»1 

three  values  has  to  be  chosen  for  h  in  order  to  make  the  order  of  (136)  a  maximum  depends 

therefore  on  (137).  h  must  be  that  one  of  the  three  numbers  p,  p+1,  p+2  for  which  the 
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X^  and  the  can  b«  chosen  to  w  to  oaks  T  a  maximum.  Let  ua  call  this  maximum 

v  u 

T^11*.  As  in  section  14  it  is  clear  that  nust  belong  to  T^11^,  as  well 

M  Tr+1'  Tr*2 . V^p-2‘ 

With  respect  to  the  two  retraining  tents  to  be  chosen  we  can  only  say  that  they  nust  be 
the  two  largest  of  the  remaining  X*s  and  tJs.  This  proves  that  the  minor  of  greatest 
order  of  nagnitude  among  all  minors  of  the  first  n-a  columns  of  &(  p)  is  in  this  case 
forned  by  the  rows 

1,2,.  .  .,pi  r+1 ,  r+2, . .  . . ,r+n-n-p-2  ,  (172) 

and  two  nore  rows,  which  must  be  those  corresponding  to  the  boundary  conditions  containing 
the  highest  order  of  differentiation  excluding  those  rows  already  contained  in  the  sequence 
(172).  A  comparison  with  part  B  of  the  third  step  of  the  Main  Theorem  shows  that  these  are 
just  the  rows  of  h[ p)  belonging  to  boundary  conditions  that  must  be  canceled  in 
application  of  that  rule.  We  know  from  the  proof  in  case  I  why  this  is  sot  the  rows  of 
A(  p)  appearing  in  the  minor  of  maximal  order  ot  the  first  n-m  columns  of  M  p)  are 
just  those  that  do  not  occur  in  the  cofactor  of  this  minor,  and  this  cofactor  determined 
the  boundary  conditions  satisfied  by  the  limit  of  U  ( x ,  p) ,  if  it  exists. 

The  exceptional  case  of  indetermination  occus,  when  the  two  additional  rows  after  the 
rows  (172)  have  been  chosen,  are  not  uniquely  determined.  We  shall  assume,  for  the 
present,  that  we  have  to  do  with  the  regular  case.  The  case  of  indetermination  will  be 
treated  in  J9, 

Our  rule  is  so  formulated  that  it  takes  also  care  of  the  case  that  X  . .  or  X  . ,  do 

pt  i  p*  . 

not  exist  because  r«p  or  r  «  p+1. 


27.  The  minor  of  maximal  order  or  magnitude  among  the  minors  of  the  first 

D 

n-m  columns  of  &(p)  can  now  be  written 


.(II) 


d£II’<0> 


H 

°  i 

V»1 


(173) 
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where 


,<II) 


H  r+n-m-H 

I  v  i 

w*1  ivrM 


(174) 


and  H  is  equal  to  p,  pH,  or  p+2 ,  as  the  case  may  be. 

As  in  case  I  we  can  then  prove  that  Q* 11  *  ?  0,  provided  assumption  9*  is  satisfied 

D 

The  cofactor  of  D* 11 ^ ( p)  will  again  be  of  the  fora 
D 


iC1*1 


where  A^11  *  is  defined  in  analogy  with  Assuaption  8*  of  the  Mein  Theorem  assures  us 

again  that 


*  0  • 


Finally,  we  find  for  A( p),  similar  to  (161) 


A(p)  -  c  b 


(id  o  i  » 

-  Vi 


in  case  IIB 


(175) 


28.  As  in  case  I,  we  could  now  calculate  the  determinants  A^{ p)  by  the  method  used  for 
Mp).  But  since  we  are  only  Interested  in  proving  that 


V  p> 

A  (p) 


«v(*.p) 


tends  to  zero,  for  v  <  n-m,  we  shall  not  calculate  the  analogue  of  (169)  for  this  case. 

As  in  proof  of  lemma  3  we  consider  first  the  A  ( p)  for  v  <  p.  For  these  v  the 

°"v 

determinant  A  (p)  does  not  contain  the  column  with  e  of  A  (p).  It  follows  from 

v  V 

lemma  1,  that  then  the  minor  of  maximal  order  among  the  minors  of  the  first  n-m-1  columns 

of  A  (p)  cannot  contain  as  factor  an  exponential  of  higher  order  than 
v 


"~T 


X 


a  7  w  -  w 

A  «  v 

• ap1 

This  is  thsrsfors  also  ths  sxponsntial  of  maximal  order  that  can  occur  in  the  asymptotic 
expression  for  A  (x).  Hence/ 


has  the  exponential  factor  e  (and  possibly  exponential  factors  of  order  1)/  and 


T7p7  °vl*'pl  '  U<p) 


has  an  exponential  factor  whose  exponent  is 


<  »<«>«« -v 


an  expression  which  tends  to  sero  as  0  ♦  «■ 

V  p> 

If  v  >  p*2,  we  prove  eimilarly  that  does  not  contain  any  exponential  factors 

(except  possibly  exponentials  of  the  order  or  magnitude  of  1 ) .  Hence  the  asymptotic 
expression  for 


yp) 

^Hv(x,p)  ,  (  v  >  pf  2 ) 


contains  an  exponential  factor  whose  exponent  is 


o  *vU)d5  /  (  v  >  p*-2)  , 


and  this  expression  tends  to  sero  as  o  ♦ 


”44' 


fcS 


29.  The  two  remaining  terse  of  (133)  require  a  more  careful  analysia.  We  treat  only  the 
case  v  -  pH,  the  caae  v  «  pt2  being  almost  identical.  Similar  as  in  section  26  we  ask 
which  choice  of  n-m-1  rows  of  A_,  .  ( p)  leads  to  the  minor  of  highest  order  among  the 

ptl 

minors  of  the  first  n-®-1  columns  of  A_^.(p).  This  minor  will  have  either 

ptl 

a\  w  a  7  w  +  w 

«  <>>  .«*  P+2 

<tf*1  IV*1 

e  or  e 

as  exponential  factor.  In  both  cases  the  asymptotic  expression  for 

A  <  p) 

AM~  Vltx'p) 

will  contain  an  exponential  factor  of  the  order  of  1.  But  the  asymptotic  expression  for 

T' 

A  (  p)  has  also  a  factor  of  the  form  a  ,  and  similarly  as  in  section  26  we  conclude 
ptl 

that  T'  is  a  sum  which  contains  the  terms  X1 ,  X^ , • • • »  X^  and  the  terms  Tf+1 , 

t  To  these  terms  one  more  term  has  to  be  added  (not  two  terms  as  in 

r+2  rtn-m-p-2 

the  case  of  A(  p) ,  because  A  ,  ( p)  must  be  expanded  in  terms  of  its  first  n-m-1,  not 

pt  1 

n-m,  columns) .  This  term  must  be  the  largest  of  the  remaining  XJs  and  tjs.  Finally, 

a  ,  i  _(H> 

A  <p)  T  '-T 

in  the  asymptotic  calculation  of  \  .  ,  we  have  to  form  g  ,  and  this  will  be 

Vi(p) 

the  power  of  o  occurring  in  the  asymptotic  expression  for  .  ,  V  J4(x,p).  This 

A  (p)  A(P)  ^ 

proves  that  the  order  of  magnitude  of  \  .  U  .  «(x,p)  is  not  greater  them  that  of 

Alp)  p+ 1 


where  s  is  the  smaller  of  the  two  terms  chosen  for  T.  7 ,  after  X..  * X* , . # • , \  and 
— — —  d  i  2  p 


ww 


,  t  _  have  been  selected,  s  may  be  any  of  the  numbers 

n-  xsr  p“  2 

V 1 '  V* 2 '  Tr+n-ra-p- 1 '  Tr+n-m-p 


In  order  to  prove  that  the  ptl  st  term  of  (133)  tends  to  zero  it  remains  only  to  show 
that  none  of  the  numbers  (176)  can  be  zero: 


*>  ■  “  X  .  In  thi»  case  the  other  on*  of  th*  two  last  numbers  chosen  for  T*11' 

P+1  b 

■uat  be  Tr+n_m.p_1>  because.  If  It  were  X^.  e  would  not  be  the  smaller  of  the  two, 

and  it  cannot  be  x  since 

r+n-*-p 

Tr+n-m-p  <  Tr+n-m-p- 1 


and  the  two  chosen  numbers  must  be  the  largest  of  the  numbers  (176).  Since  s  must  be  the 
smaller  of  these  two  numbers,  we  conclude  that 


On  the  other  hand. 


X  <  x  .  , 

pH  r+n-m-p-1 


V+1  >  Tr+n-m-p 


I<  X  “  0.  X  ,  must  be  the  last  of  the  numbers  X,  hence 
P+ *  p+ ■  1 

p+1  “  r  . 

It  then  follows  from  (177)  that  x  cannot  exist,  in  other  words 

r+n-m-p 

x  .  .  «  x 

r+n-m-p-1  n 


(178)  and  (179)  imply 


r+n-m-p-1  -  n 


m  ■  0  , 


a  case  excluded  from  our  considerations.  Hence,  X  .  .  >  0. 

p+1 

b)  a  «  X^,.  Similarly  as  in  a)  it  follows  that 
P+2 

^p+2  >  Tr+n-m-p-1 

X  . ,  ”  0  would  imply  p+2  -  r  and  r+n-m-p- 2  »  n,  hence 
P+* 

ra  ■  0  . 


c)  and  d)  a  -  x  ,  or  s  "  x  .  .We  leave  the  proofs  in  these  cases  to 

rTn“m*p'*  1  r+n*B“p 

the  reader,  since  they  follow  exactly  the  pattern  of  the  proofs  in  a)  and  b)  with  the  same 
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30.  Thus  we  have  finished  the  proof  of  the  following  theorem. 
Theorem  3i  If  the  assumptions  1*  -  9*  of  section  3  are  satisfied,  if 

0  <  r-p  <  mv  2  , 


and  if 


n-m  s  0  (nod  4),  and  b^  <  0 
or 

n-m  =  2  (mod  4),  and  bQ  >  0  , 

then,  as  p  ♦  «,  the  solution  U(x, p)  of  the  problem  defined  by  (101)  and  (108)  converges 
in  a  <  x  <  8  to  that  solution  of  the  limiting  differential  equation  M(y)  -  0  which 
satisfies  the  m  boundary  conditions  not  canceled  in  application  of  the  rule  of  the  Main 
Theorem  for  the  case  (xv),  unless  we  have  to  do  with  the  case  of  indetermination. 


31.  The  case  IIC.  In  case  XIC  (i.e.,  if  r-p  >  mt2)  one  can  again,  as  in  section  23, 
consider  the  (r-m)th  term  of  the  right  lumber  of  (133)  and  prove  that  its  asymptotic 
expression  contains  the  exponential  factor 

<>{/%<  Od5  -  w  ) 

a  r-w  r-m 

0  • 

This  is  sufficient  to  show  that  this  term  must  tend  to  infinity,  since  this  exponential 
tends  to  infinity  more  strongly  than  any  potter  of  a  may  tend  to  zero. 


§9.  The  Case  of  Indetermlnatlon. 

32.  let  us  assume  that  we  cannot  decide  in  a  uniquely  determined  way,  which  is  the  last 
row  of  A(  p)  to  be  chosen  in  order  to  obtain  the  minor  of  greatest  order  of  magnitude 
among  the  minors  of  the  first  n-m  columns  of  &( p) .  This  is  the  case  which  we  have  left 
aside  in  section  26.  We  have  seen  there  that  this  occurrence  means  that  we  have  to  do  with 
the  case  of  indetermination  of  the  Main  Theorem. 

In  this  case  we  must  have 

Xt“  \ 

where  i  is  one  of  the  numbers  pM,  p+2  and  k  one  of  the  numbers  r+n-m-p- 1 ,  r+n-m-p. 
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Furthermore,  X^  and  must  ba  tha  aacond  and  third  in  aisa  of  the  nunbera  and 

■t.,  excluding  tha  numbers  X.,X_,...,X  and  t  .,  _.  For  simplicity 

i  id  p  r+i  r+2  zrn*Bhp"2 

lat  ua  assume  that  t  *  pM.  {Tha  raaaoning  ia  the  same  for  t  -  pf  2. )  Denote  by  D1  ( p) 
and  ( 0)  tha  ainora  of  the  firat  n-ra  columns  of  tha  determinant  At p) ,  with  tha  rows 

1.2..  .  .,p,  p+1»  r+1,  r+2,...,k-l 

1.2. . ..,p;  r+1,  r4-2,...,k-1,  k  , 

raapactivaly.  Than  D^tp)  and  ( p)  are  both  minora  of  the  same  maximal  order,  and  we 
cannot  raaaon  aa  in  aection  16.  But,  uaing  formula  (138)  we  can  write 


ow 


D^p)  *  F(  p)  e 


P*1 


[d,] 


D2(p)  -  Tip) 


td2)  , 


where 

^  *v  +  I  t  0!  w 

Ft  p)  -  O'"1  ^  “a  -1 

d1  and  d2  are  certain  conatanta,  which  can  eaaily  be  proved  to  be  different  from  zero  by 
the  method  of  aection  IS,  provided  assumption  9*  is  satisfied  for  each  of  the  two  ways  of 
canceling  corresponding  to  D^tp)  and  to  02(p),  respectively. 

The  cofactors  of  the  minors  D ^  { p)  and  D2<p)  in  At  p)  can  be  written  in  the  form 
[5j]  and  [Aj],  where  6^  and  A2  are  certain  constants.  If  assumption  8*  is 
satisfied  for  both  waya  of  applying  the  cancellation  rule  of  the  Main  Theorem,  then  we  can 
be  sure  that 

A,  f  0,  A2  +  0  . 

expanding  At  p)  in  terms  of  its  first  n-m  columns  we  find  then 


At  p)  -D^pMA,!  tD2tp)lA2]  -  Ft  p)  {e  d^]  ±  ld2A2l}  .  (180) 
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33.  While  in  all  cases  so  far  the  last  a  terms  of  the  sum  (133)  were  convergent, 
the  convergence  or  non-convergence  of  U(x, p)  depended  entirely  upon  the  first  n- 
of  (133),  in  this  case  the  last  m  terms  of  (133)  will,  in  general,  be  divergent. 
In  fact,  let  n-m  <  j  <  n.  Then  we  find  immediately  in  analogy  with  (100), 

oh 

4j(P)  *F(p){e  ^’td^j)  +  (d2«2jn  ' 

where  [6,^]  and  (5^1  are  the  determinants  obtained  instead  of  (6^)  and  S2), 
j-th  column  of  A(  p)  is  replaced  by 


t 


(180),  (181)  and  (115)  show  that  the  j-th  (J  n-m)  term  of  (133)  is  equal  to 

ow 


Mp> 


d(rt  vx'p) 


>1 


Of 


lVlj5  +  [V21>  .  .  .  . 

- * - J-  [u  (x)  ] 


P*'1  [d,  6,]  +  [d2S2)  3"n"B 


and  that 


n  A.  (  p) 


OW  n  n 

«  l  [dt«  u  (x)  )  +  l  [A2&2  u  (x)  1 

i**n-m _ 3  j-n»ra  j-n+a _ 3  j-n-m 


0*»  . 

e  ^  Id,  6,1 


Ed252l 


Because  of  the  oscillating  factor  e  pf1  this  expression  can  converge  only  if  the 
determinant 


i  d,a,ju  <*> 

j-n-m  J  j-n-o 


Vi 


l  d2«2  u  (x) 
j-n-m  3  j*n+m 


V2 


vanishes,  i.e.  if 


so  that 
m  terms 

(181) 

if  the 


(102) 

(183) 
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n  42j| 

l  u  (x)  -  0  . 

j-n-m  «1  i2  [  j-n+m 

Sine*  the  uy(x),  (u  -  1.2, ...,m)  are,  by  assumption,  linearly  independent,  thie 
leads  to  the  n  conditions 

$25ij  “  5^2j  “  °  ^  •  (184) 

The  left  sides  of  (184)  depend  on  the  numbers  t,,  ,  (184)  represents  therefore  a 

12  n 

set  of  conditions  on  the  prescribed  boundary  values e  As  in  the  previous  cases  we  assune 
that  the  tjs  do  not  have  the  very  special  values  required  by  (184).  Then,  the  value  of 
(183)  oscillates,  as  a  *  •,  but  remains  bounded. 

As  to  the  first  n-m  terms  of  the  sum  in  (133),  it  can  be  proved  exactly  as  in  the 
regular  case  IIB  (section  27-29)  that  they  all  tend  to  zero. 

This  completes  the  proof  of  the  non-converqence  of  U(x,p)  in  the  case  of 
indetermination  and  also  the  proof  of  the  whole  Main  Theorem. 


34.  Remark i  Going  over  our  whole  proof  we  see  that  the  assumption  m  >  0  was  not  used  at 


all  in  the  proof  for  case  I,  i.e.,  for  the  cases 


©-  © 


of  (113).  Hence,  all  our 


In  case  II,  i.e.  (IVJ,  t*1®  hypothesis  m  >  0  was  used  only  to  exclude  a  rather 


results  in  these  cases  remain  valid  for  m  ■  0  also. 


special  occurrence  in  section  29.  it  would  not  be  difficult  to  formulate  a  general  theorem 


for  the  case  m  -  0  also. 
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FURTHER  RESULTS  IN  THE  CASE  n-m  -  1 
$1 .  Generalization  of  the  Boundary  Conditions. 

35.  In  the  special  case  n-m  -  1  It  is  not  difficult  to  replace  the  boundary  conditions 
(108)  by  the  more  general  ones 


n  (v-1)  n  <  w-1  > 

I  oivy  (a)  ♦  I  8  y  (8) 

v-1  u-1 


(1  m  1/2/ .. .  ,n) 


(201) 


where  the  constants  ci^,  8^  and  are  only  restricted  by  the  condition  that  the  n 

boundary  conditions  are  independent  and  compatible. 

We  note  first:  If,  in  (201),  the  L, (y.  _ .  .. _  .  _  _  .  .  .  .  .  .. 

i  J)  and  the  are  subjected  to  the  same 

linear  transformation  with  constant  coefficients  and  non-vaniBhing  determinant,  then  the 
resulting  equations 


n  ( v-1 )  n  ( vi- 1 ) 

LMy)  -  l  a'  y  (a)  +  J  B*  y  (0)  -  « j  (202) 

v-1  v-1  W 

constitute  a  set  of  boundary  conditions  equivalent  to  (201)  in  the  sense  that  a  function 
satisfying  (201)  satisfies  also  (202)  and  vice  versa. 

From  this  remark  we  see  that  we  can  assume  without  loss  of  generality  that  not  all  the 
or  all  the  are  zero,  because  in  that  case  the  boundary  conditions  (201)  would 

be  equivalent  to  ordinary  initial  conditions,  which  are  a  special  case  of  boundary 
conditions  of  the  type  (108). 

36.  Denote  by  s  the  greatest  value  of  v  for  which  at  least  one  of  the  a,  is  not 

lv 

zero,  and  by  t  the  greatest  value  of  u  for  which  at  least  one  of  the  B^  is  not  zero. 
There  exists  a  non-degenerate  linear  transformation  with  constant  coefficients  transforming 
(201)  into  the  system  of  equivalent  boundary  conditions 
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(203) 


o*i 


n 

■  I 

vl 


a. 

a  i  v 


<  v-1> 
y  (a) 


l 

ifl 


a  ip 


C  M-1 » 
y  ( 8)  - 


A 


having  th«  property  that 

a.  i*  0.  cu^  ■••••*  a  “0  . 

a  It  r  ait  a  ne 

Such  a  transformation  can  be  chosen  in  many  ways.  Similarly,  there  is  a  transformation 
changing  (201)  into 


6 


n  (\>-1)  n  <p-1) 

*1  S  J,  0°ivy  (a)  +  J,  06iyy  (6> 


)*■  1 


0li 


(204) 


having  the  property  that 

S61t  *  °'  002t  BKnt 

For  n-m  *•  1  aquation  (114)  and  (IIS)  reduce  to 


0  . 


-p  £y°d* 

U^x.p)  -  e  (n(x)l 


U  (x.p)  -  [uU>,  ,  (p-  1,2,...,n-1)  . 


Let  us  assume  first  that  b0(x)  >  0.  Then  we  use  (203)  instead  of  (201)  and 


immediately  that 


yiy  -  ±  p*" 1  ( bQ ( o0 1  •  n(o>] 


is  of  greater  order  of  magnitude  than  all  the  other  L^(U^) ,  (i  “  2,3,...,n).  Solving 
equations  (125)  asymptotically  in  this  case,  by  the  method  used  in  the  first  chapter,  we 
see  that  we  can  reason  exactly  as  there.  Of  the  assumptions  of  the  Main  Theorem  we  need 
only  1*,  2*,  3*,  5*,  6*  and  an  assumption  corresponding  to  8*  which  states  that 


aVu1>' 


'  _L2<“n> 
a  2  n 


L  (u  ), 
an  1 


*  •  > 
an  n 


+  0 


(205) 


We  see  then  easily  that  the  solution  of  (101)  and  (201)  tends,  for  n-m  *  1,  and  bp  >  0, 

to  that  solution  of  M(y)  ■  0  which  satisfies  the  boundary  conditions 

L, (y)  -  t.  ,  (1  -  2,3,. ..,n)  (206) 

a  1  a  1 


but  not  the  boundary  condition  L  ( y )  -  1.  except  for  special  values  of  the 

a  l  ex  i 

a  similar  reasoning  in  the  case  bg(x)  <  0  we  find  that  in  that  case  u(x)  ■ 
satisfies  M(y)  -  0  and  the  boundary  conditions 


i’s.  Using 
lira  U(x,  p) 

P-m. 


0Vy)  ■  8*i  •  (i  “  2'3 


(207) 


provided 


V 


BL2(u1)'  *  *  *  '  bW 


eLn(U1)'  •  •  •  '  bW 


+  0 


(208) 


This  result  can  be  formulated  in  a  somewhat  more  symmetrical  form.  TO  that  end  note 
that  the  boundary  conditions  (206)  do  not  involve  any  more  the  highest  derivative  at  x  ■  a 
occurring  in  (201).  any  linear  combination  of  the  equations  (201)  which  does  not  contain 
(s-1) 

y  (a)  must  be  linearly  dependent  on  equations  (206)  and  no  linear  combination  of 

(s-1)  (t-1) 

equations  (206)  contains  y  (a)  •  Similarly  for  (207)  with  respect  to  y  ( 8)  •  Hence,  we 

can  state  the  following  theorem. 


-53- 


Theorem  4»  if  n-m  “  1  and  if  tha  conditions  1»,  2*,  3*,  5*,  6*  of  tha  Ha  in 
Thaorea  of  chaptar  I  as  wall  aa  (205)  and  (208)  ara  aatiaf lad,  than  tha  function 
U(x, p)  satisfying  (101)  and  tha  boundary  conditions  (201)  tanda  with  incraaaing 
p  to  a  solution  u(x)  of  *K  y)  “  0.  According  as  bg  >  0  or  b^  <  0 
tha  function  u(x)  satisfies  all  boundary  conditions  that  dapand  linearly  on  (201) 
and  do  not  contain  tha  highest  derivative  at  x  -  a  or  x  “  0,  respectively, 
occurring  in  (201). 

Remark;  The  conditions  (205)  and  (208)  can  be  formulated  in  a  way  independent  of  tha 
particular  choice  of  the  fundamental  systaa  u^(x),  ( u  •  1,2,...,n-1) ,  by  saying  that  we 
assume  that  only  the  function  u(x)  =0  satisfies  tha  differential  equation  M(y)  «  0  and 
the  homogeneous  boundary  conditions  corresponding  to  (206)  or  (207),  respectively. 

Example :  n  -  3,  a  -  2, 

y(  a)  -  y'  (  a)  +  y*  •  (  8)  -  i, 

y(n)  +  y*  (a)  +  y(0)  -  y*  ( 0)  -  *2 

y«(a)  -  2y'  ’  (  8)  -  t3 

If  bQ  >  0,  then  lira  U(x,  p)  satisfies  the  boundary  conditions 
P*“ 

2y(  a)  +  y(  0)  -  y*  (  8)  ♦  y*  ’  (  0)  -  i,  +  t2 

y(a)  +  y(0)  -  y’{0)  +  2y"(8>  -  »2  "  »3  • 

But  if  bQ  <  0,  then  lim  U(x,p)  satisfies  the  boundary  conditions 

p+- 

y(  a)  ♦  y*  ( a)  ♦  y(  8)  -  y’  (  8)  -  *2 
2y(a>  -  y'(a)  -  2t,  +  t3 
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$2.  The  "Stretching*  of  the  Boundary  Um 
37.  For  the  relatively  simple  types  of  boundary  layer  problems  with  which  this 
investigation  is  concerned  we  have  been  able  to  develop  a  method  that  allows  us  to 
calculate  asymptotic  expressions  for  the  solution  of  the  boundary  value  problem  (compare, 
e.g.,  formula  (171)).  From  these  asymptotic  expressions  one  can  easily  obtain  all  desired 
information  about  the  behavior  of  the  solution  of  the  boundary  value  problem  near  the 
endpoints  for  large  values  of  p.  (Compare,  e.g.,  the  corollary  in  section  20). 

In  the  more  complicated  boundary  layer  problems  occurring  in  physics  such  complete 
asymptotic  solutions  are  often  not  available.  In  those  cases  it  is  customary  to  transform 
the  given  boundary  value  problem,  by  a  change  of  the  independent  variable,  into  a  new 

boundary  value  problem  which  does  not  tend  to  a  problem  of  lower  order  when  p  tends  to 

infinity. 

as  an  example  for  such  a  transformation  we  take  the  differential  equation  (101),  for 
the  special  case 

bg(x)  >  0  . 

The  case  bQ(x)  <  0  can  be  treated  analogously.  Without  loss  of  generality  we  may  further 
assume  that 

a  -  0  .  (209) 

We  shall  refer  to  this  boundary  value  problem  as  the  problem  (L) . 

We  now  introduce  the  new  independent  variable 

s  *  (Sc  ,  (210) 

and  transform  the  boundary  value  problem  (L)  into  an  equivalent  problem  in  x,  to  which  we 
shall  refer  as  the  problem  <L).  Let  U(x,p)  be  the  solution  of  the  problem  (L) .  0(x,p) 
or  some  of  its  derivatives  will  have  a  boundary  layer  at  x  «  a  »  0.  The  function  U(x, p) 
is  changed,  by  the  transformation  (210)  into 

U(x,  p)  -  U(£,  p)  .  (211) 

P 

U(z,p)  is  the  solution  of  the  problem  (iT).  Since 

U<i)(x,p)  -  p_i  U(1)  (x,p)  ,  (212) 
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the  problem  (L)  can  be  written 


fir)  5  y(B>  ♦  b0(2)y(*)  ♦  l  ev(i)p-Vn-V>  ♦  l  **» 

p  v*l  p 


W  P 


(2.13) 


(X.)  -X. 

y  1  (pg)  -  tlP 


0  <  1  <  r 


L1  S 


<  V 

y  (0)  -  ttp 


(214) 

r+1  <  i  <  n  . 

If  we  let  p  tend  to  Infinity  In  the  coefficients  of  (213)  we  obtain  the  simple 
"limiting"  differential  equation 

y(n>  +  b„<0)  y(B)  -  0  .  (215) 

It  may  be  expected  that  the  function  U(x,  p)  tends  with  increasing  p  to  a  function 

u(s)  -  11m  <7(s,p)  (216) 

p-M 


which  satisfies  the  differential  equation  (215).  Since  (215)  is  of  the  same  order  as  (213) 
we  expect  that  the  function  V(x, p)  will  not  have  a  boundary  layer  for  large  p. 

The  transformation  (210)  may  be  described  as  a  stretching  of  the  function  U(x,p).  If 
U(s,p)  does  not  have  a  boundary  layer,  we  have,  in  a  way,  "stretched  out"  the  boundary 
layer . 


The  problem  arises  than  what  boundary  conditions  are  satisfied  by  the  limiting 
function  u(z).  If  the  boundary  layer  of  the  function  U(x,p)  at  x  ■  0  occurs  in 
U(x,p)  itself  and  not  in  a  derivative  of  U(x, p),  the  interpretation  of  (210)  as  a 
stretching  which  becomes  infinite  when  p  ♦  «•,  suggests  that  we  have 

u(«)  -  u<0)  , 

where  u(0)  is  the  value  assumed  by  the  function 

u(x)  -  11m  U(x,p)  (217) 

p*m 

at  x  «  0.  u(0)  will,  in  general,  be  different  from  the  boundary  value  prescribed  for 

U(x,p)  at  x  •  0. 
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Tha  transformation  (210)  is  frequently  used  for  tha  solution  of  mors  complicated 
boundary  layer  problems.  Tha  following  points  are  than  usually  taken  for  granted  without 
proof: 

(a)  That  lim  U(z,p)  exists. 

P-*** 

(b)  That  the  limit  u(z)  satisfies  the  limiting  differential  equation 

(c)  That  u(«)  -  u(0). 

In  our  investigation  we  have  been  able  to  find  an  asymptotic  approximation  for 
0(x,p)  directly,  so  that  we  did  not  need  the  transformation  (210).  But  we  are  now  able 
to  prove  the  statements  (a),(b),  (c),  for  our  problem  (L).  This  is  what  we  are  going  to  do 
in  this  $. 


38.  In  chapter  I  we  have  derived  for  U(x,p)  the  following  asymptotic  representation 
(compare  (171)): 

'V)  P  fo 

U(x.p)  -  (sn(X)Jp  e  ♦  (u(x)]  .  (218) 

In  this  formula  we  are  using  the  following  abbreviations: 


r  w  m  \  \  Tr+1  *1 

*  ’  f-b0<0)>  T 


(219) 


where 


i  -  (-D 


r+1 


vv  • 

•  •  •  •  L,(u  ) 

1  B 

VV  • 

....  Lr(Un) 

WV 

•  •  •  •  wv 

w  • 

(220) 


( Compare  formula  (159)), 
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•  •  •  » 


*1  !.,<«,)  ....  L,^) 

*2  W  *  *  *  *  W 


t  I.  (U,)  .  .  .  .  L  (tt  ) 
n  n  <  n  ■ 


(221) 


(compare  section  18); 


and 


*(x)  -  -bQ{x) 


(222) 


u(x) 


m  5, 

V  _!±1! 
L.  6 
irl 


«p(x) 


(223) 


(Compare  section  17  for  the  6^^,  -  we  have  dropped  the  Index  b,  used  there,  as 

unnecessary  here  -/  and  compare  section  6,  theorem  1  for  the  u^{x)). 

In  consequence  of  theorem  1,  formulas  (118)  and  (120),  we  may  differentiate  (218) 
formally  at  least  n-t  times,  l.a. 


0(1)(x,p) 


1-T 


•  [*  (x)n(x))p 


r+i 


£  *<0*5 


[uU)(x)J, 


(1-0, 1,...,n-1) 


(224) 


39.  Knowing  0(x,p)  and  Its  derivatives  ws  can  now  easily  calculate  the  function  U(r, p) 
and  its  derivatives  with  respect  to  s.  For,  considering  that 

a  £/p  *(C)d5  /?  *<-£)di; 

a  0  -  e  0  (225) 

we  find,  upon  substitution  of  (212)  into  (224),  that 


U,l)(s.p) 


tiW>l(f)  «n(*;)I  P 
P  P 


r+1 


f! 


p~l<£)],  (1  -  0, 1, . . .  ,n-1) 
P 


(226) 


Thus  we  have  solved  the  problem  (L ) . 
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The  result  of  the  passage  to  tha  Halt  In  (226)  can  bs  moat  easily  axprsssad  by  ons 


formula  comprising  all  cases  that  can  arise.  If  we  Introduce  the  following  symbol i 

(  0.  If  t  ^  0 

e(t)  -  {  (227) 

l  1,  if  t  -  0  . 

Then  we  can  write 

u(s)  -  11a  0(a,  p)  -  (i*W(0))e2^<0)e(Tr(,1)  ♦  u(0)  .  (228) 

P*» 

Formula  (228)  is  valid  for  0  <  s  <  «. 

We  can  now  confirm  the  three  unproved  statements  (a),  (b),  (c)  of  section  37.  u(s) 
does  exist,  and  it  satisfies  the  limiting  differential  equation  (215),  as  we  may  readily 
verify  by  substitution.  In  fact,  u(s)  reduces  to  a  constant,  unless  1  *  0.  Finally 
we  see  from  (228)  that  u(  -)  -  u(0) .  This  proves  the  stateaents  (a),  (b),  (c). 

For  the  limit  of  the  i-th  derivative  of  U(x,p)  we  find,  from  (226), 

limO<i,(*.p)  -  U*l(0))e*'',<0>c(TrM)  ,  (1  •  1,2,...,n-1)  ,  (229) 

p+m  1 

which  is  also  the  i-th  derivative  of  u(s).  We  conclude  from  (228)  and  (229)  that  0(0, p) 

J  4  \ 

and  0  (0,p)  remain  finite  as  p  ♦  in  other  words,  0(*,p)  does  not  have  a  boundary 

layer  at  x  »  0,  the  boundary  layer  has  been  'stretched  out*. 

40.  We  have  seen  that  the  limiting  differential  equation  (215)  is  satisfied  by  the 
function  u(x)  obtained  by  passing  to  the  limit  in  the  solution  U(z,p)  of  the  problem 
(L) .  How  can  we  find  a  complete  set  of  n  boundary  conditions  satisfied  by  u( s) 7 

If  we  formally  let  p  tend  to  infinity  in  the  boundary  conditions  (214),  we  obtain 
boundary  conditions 


(->  -  e(Xi) 


(i  -  1,2 


r) 


(230) 


y 


(0)  -  e(Tt) 


(1  -  r+1,  r+2, . . . ,n) 


(231) 


(Ti) 


From  the  remark  made  at  the  end  of  the  laat  section  it  follows  immediately  that  u(z) 
satisfies  the  boundary  conditions  (231).  If  »  a,  the  last  boundary  conditions  (230) 
can  only  be  satisfied  if,  by  coincidence,  i »  u(0).  We  therefore  replace  the  last 
boundary  condition  of  (230)  by  the  condition 

u(»)  -  u(0)  , 

which  we  have  proved  to  be  satisfied.  The  other  boundary  conditions  (230)  are  certainly 
satisfied,  since  all  the  derivatives  of  u(s)  vanish  at  z  •  *. 

$3.  The  Won-Homogeneous  Differential  Equation  —  S(y)  ♦  M(y)  **  f(x). 

41.  Introductions  It  is  an  open  question  whether  the  Main  Theorem  remains  valid  in  full 

generality  for  the  non-homogeneous  differential  equation 

-  H(y)  +  M(y)  -  f(x)  .  (232) 

P 

But  we  will  be  able  to  answer  this  question  in  the  affirmative  when 

n  -  m  -  1  ,  (233) 

provided  none  of  the  boundary  conditions  not  canceled  in  application  of  the  rule  of  the 
Main  Theorem  involves  an  order  of  differentiation  greater  than  m-1.  The  meaning  of  this 
latter  condition  is  easily  understandable i  If  one  of  the  uncanceled  boundary  conditions  is 
of  the  order  n-1,  then  the  boundary  value  problem  formed  by  the  limiting  differential 
equation 

M(y)  -  f(x) 

and  the  uncanceled  boundary  conditions  is  of  a  type  to  which  the  usual  method  of  solution 
by  means  of  the  Green's  Function  cannot  be  applied,  since  this  method  presupposes  that  the 
boundary  conditions  are  of  lower  order  of  differentiation  than  the  differential  equation. 

In  our  treatment  of  the  homogeneous  differential  equation  the  relative  order  of 
differentiation  of  the  limiting  differential  equation  and  of  the  remaining  boundary 


-60 


conditions  did  not  play  any  important  role,  and  it  is  by  no  moans  certain  that  the 
condition  above  is  really  necessary  in  the  non-homogeneous  case.  But  it  simplifies  our 
proof  greatly. 

42.  The  problem  defined  by  (232)  and  (233)  and  the  boundary  conditions  (108)  will  be 
called  the  problem  (N).  If,  instead  of  (108),  we  prescribe  the  corresponding  homogeneous 
boundary  conditions,  we  shall  speak  of  the  problem  (N ' ) .  We  assume  that  0  <  r  <  n,  i.e. 
we  consider  only  actual  boundary  conditions  leaving  aside  the  Initial  value  problem.  (The 
initial  value  problem  can  be  treated  by  the  same  method.)  It  is  easy  to  extend  the  proof 
below  to  the  more  general  boundary  conditions  (201). 

Let  z( x , p)  be  the  solution  of  (W),  *(x,p)  the  solution  of  (N ' )  and  U(x,p>  the 
solution  of  the  homogeneous  differential  equation  (101)  satisfying  the  non-homogeneous 
boundary  conditions  (108).  Then 

Z( x , p)  -  z(x,p)  +  U(x,p)  .  (234) 

Since  the  asymptotic  behavior  of  U(x,p)  for  large  has  already  been  investigated  in 
chapter  I  it  is  sufficient  to  discuss  the  problem  (N'). 

43.  The  Green's  Function: 

It  is  known  that  the  function  z(x,p)  can  be  written  in  the  form 

z(x,p)  »  f®  G( x, t ,  p)  p  f(t)  dt  .  (235) 

Cl 

The  "Green's  Function"  G(x,t,p)  can  be  constructed  in  the  following  manner: 

Let  y,(x,p),  y„(x,p)  •  •  •  y  (x,p)  be  a  fundamental  system  of  solutions  of 
*  2  n 

—  N(y)  +  M(y)  «  0,  and  set 
P 
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Mt,p) 


y(n-1)(t  p)  y(n  .  ,  .  y^n  ^(t.p) 

y‘n_2)(t,p)  y<n'2)(t,p)  •  •  .y^n_2)(t,p) 


y^t.p)  y2(t,p) 


yn(t.p) 


k(x,t,p)  -  agn(x-t) 


y^x.p) 


y2(*»p> 


•  •  •  y_<x»p> 

n 


yin  2)(t,p)  yin  2)(t,p)  •  •  •  y*n  2)<t,p) 


y^t/p)  y2(t,p)  •  •  •  yn(t,p) 


(236) 


and 


g(x,t,p) 


1  k(x,t,p) 

2  h(t,p) 


Than 


G(x,t,p) 


(-1) 


H(x,t,p) 
A<  p) 


(237) 


(238) 


whare 


and 


H(x,t,p) 


y^XfP)  y2(x»P> 

W  W 


y  (xrp)  g(x,t,p) 
n 

W  (,  (g) 


(239) 


w 


L  (y,) 
n  2 


\  (y„> 

n  n 


L  (g> 
n 
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J 


Vy1>  Vy2J  •  •  *  W 


l2(y,>  i.2(y2)  •  ■  •  L2(yn) 


l  <y,)  i*  (y  )  •  •  •  l  <y  ) 
n  i  n  2  n  n 


i  “  1,2, ...,n,  means  that  the  operator  4  ia  applied  to  g(x,t,p)  considered  a 


function  of  x. 


44.  The  asymptotic  value  of  the  Green's  Function: 


We  take  as  the  fundamental  system  y^fx, p)  the  functions  (114),  (115)  specialized  for 


n-m  ■  1 ,  i .  e . 


p  r  v>(odc 

y^x.p)  -  v(x,p)  -  e  (n(x)] 


y1+u<x,P)  "  V*,p)  ”  luy<*)! 


(p  -  1.2 . n-1 ) , 


*<x)  -  -  b0(x)  . 

(Note  that  our  notations  differ  somewhat  from  those  of  chapter  I  and  of  chapter  II,  §1.) 


We  have  then 


L[(V), 


1  Li(k> 

V9)  "  2  h(t,p) 


4(0,). 


Ll<Vi> 


(n-2)  (  n-2 )  „(n_2)(a.  . 

L  ( k)  -  V  (t<p)  °1  (t'p)  *  •  *  Vi  {t'p) 


V(t,p) 


o, (t, p)  ...  a (t,p) 

I  n 


If' ^  f’  >•  >  . 


4 


with 


Tjet  us  assume  that 


Li  -  L1  • 

for 

i  <  r 

L*--Li  , 

for 

1  >  r 

h0(x) 

<  0 

* 

Then 


L’(V) 


X.  X. 

P  tv>  1(  B)  ti<6))epw  ,  for  i  <  r 


Ti  Ti 

.  -  p  If  ( a)  n<  a)  ]  ,  for  1  >  r 


where 


w  -  rB  *<£)<*£ 

Cl 


If  we  expand  (242)  in  terms  of  its  first  column,  we  see  that  in  a  <  t  <  8  the 
etc.  terms  of  the  expansion  are  of  lower  order  than  the  second  term.  Hence 

1^00  -  t-'(V)[h(t)l  -  V(n_2)(t,p)  tl^OO)  , 


where 


h(t) 


u<n-2»(t)  u<"'2>(t)  •  •  •  u(n'2)(t) 

I  i  n~*  l 

u<n"3)(t)  u‘n'3)(t)  •  •  •  u(n‘3)(t) 

I  2  n- 1 


U?(t) 


u2(t) 


•  •  U  (t> 

n-1 


and 


k(x,t)  -  sgn(x-t) 


ut(x) 


uin-3,(t) 


ll  *  •  I  U  (  V  ) 

2  n-11  ' 


(n-3 ),  n-3 

J,  (t>  •  •  •  u  ,  (t) 

2  n-1 


u^t) 


u2(t) 


•  *  ‘  Vl(t) 


(243) 

3rd,  4th, 

(244) 

(245) 

(246) 
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Expanding  h(t,p)  in  tarna  of  its  first  column  we  obtain  similarly 

h(t,p)  -  Vtn-1)(t,p)  (h(t)],  for  a  <  t  <  0  .  (2 

h(t)  i1  0,  because  otherwise  ^(t),  Uj( t) , . . . ,un_j( t)  would  be  linearly  dependent. 
From  (237),  (243),  (244)  and  (247)  we  find 


Lt(g(x,t,p) ) 


2Vln-1)(t,p) 


g(x,t) 


X(x,t) 


2h(t> 

Furthermore,  expansion  of  (236)  in  terms  of  its  first  column  yields 

k(x,t,p)  -  sgn(x-t)  V(x,p)  [h(t))  -  V<n-2)(t,p)  [k(t>] 
and  therefore,  because  of  (237)  and  (247) 


g(x,t.p)  -  sgn(x-t)  — - —777 

2V(n1)(t,p)  p  {t) 


(g(Xft)i  . 


Now  we  substitute  (250)  and  (248)  into  (239)  and  see  that 

H(x,t, p)  -  H 1  (x,t, p)  +  H2(x,t,p) 


H  (x,t,p)  -  - - - 

’  2V(n-1,(t,p) 


V(x , p) 

(u ^ ( x ) ]  »  •  •  tu|l_1(x)l,  sgn(x-t) 

V(x,p) 

X.,(V) 

•  •  • 

(L1(u1>)  •  •  •  [Lt(Un_t)] 

x.;(v) 

•  • 

•  •  • 

I.  (V) 
n 

[t  (U  )I  •  *  ■  [L  (U  )] 

n  i  n  n- » 

•  • 

I»'  (V) 
n 
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V(X,P) 

lu^x)] 

•  •  •  lu  ,  ( X )  ] 
ft—  1 

Ig(x,t)) 

I»1  (V) 

•  •  • 

Cl*1<u1 )  ] 

•  •  •  [L  (u  )1 

1  n—  i 

tx.1(g)l 

•  •  • 

(253) 

•  •  • 

t  (V) 

ft 

•  •  •  •  • 

(VVi 

lWi» 

s  s  • 

(Ln(i)] 

The  asymptotic  expression  for  A(p)  Is  (compare  chapter  X) 

X. 

A(p)  -  L1(V>[A]  -  (*,($>  n(S)  AJepw 


W  •  *  *  ^2(un-1 ^ 


L  (u  )  •  •  •  L  (u  ,) 
n  I  n  n- l 


As  In  chapter  I  we  make  the  assumption 

A  *  0  .  (256) 

(A  is  identical  with  ^  of  formula  (159)  for  this  particular  case.) 

45.  Corresponding  to  the  representation  of  H(x,t,p)  as  a  sum  of  two  terms  in  (251)  we 
find,  upon  substitution  of  (251)  into  (238), 

G(x,tp)  -  G^(x,t,p)  +  G2(x,t,p)  (257) 

with 


Gj(x,t,p)  -  (-1) 


n  Hj(x,tp) 


H  (x,t,p) 

n  i 


G2(x,t,p)  -  (-1) 


Using  (235)  this  leads  to  a  representation  of  the  solution  s(x,p)  as  a  sum  of  two 
integrals: 


--  ■*  ■'  •''-‘i.  *  .» 
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aw**-* 


(260) 


with 


*(x,p)  -  z  ^  (x ,  p)  «•  z2<x,p) 


z^x.p)  -  /®  G1(x,tp)  p  f  ( t )  dt  (259) 

z,(x,p)  -  JB  G  (x, t ,  p)  p  f(t)dt  .  (260) 

4  Cl  2 

We  shall  prove  that,  as  p  tends  to  infinity,  z^(x,p)  tends  to  zero,  while  z2(x,p) 
tends  to  a  solution  of  M(y)  -  f(x). 


46.  In  this  section  we  are  going  to  show  that 

lira  z .  ( x,  p)  -  lira  fBG,(x,t,p)  p  f(t)dt  -  0  .  (261) 

p-M»  p-M» 


To  this  end  we  write  (259)  in  the  form 


z1(x,p)  »  /*  G1 (x,t, p)  p  f(t)dt  ♦  /B  G^ (x,t, p)  p  f(t)dt 


(262) 


and  prove  that  each  of  the  two  integrals  in  (262)  tends  to  zero. 

a)  In  the  first  integral  of  (262)  we  have  t  <  x.  Hence,  G^(x,t,p)  has  to  be 
determined  with  +  V(x,p)  as  the  last  terra  of  the  first  row  of  the  determinant  in  (252). 

In  order  to  find  an  asymptotic  expression  for  H^(x,t,p)  we  expand  the  determinant  in 
(252)  in  terras  of  the  minors  formed  by  its  first  and  last  columns.  These  minors  are  either 
zero  or  of  the  form 


±  2  iyv)  Lg(V> 


i  2  V(x ,  p)I«g(V) 


T  <  r 
e  >  r 


i.e.  they  are  of  the  form 


\  +T 

.  T  e 


tql 


T 

e 

P 


Ui 

o  |*  f(5)dC 


fa'  I 


(263) 


(264) 
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are 


where  q  and  q*  are  certain  numbers  different  from  zero.  Since  and  t^+  ^ 

greater  than  all  the  other  X|s  and  t|s,  it  follows  from  (264)  that,  in  a  <  x  <  B,  the 
expression  of  highest  order  of  magnitude  among  all  the  minor  (263)  is  ±  2  Ll(V)  L^V). 
and  no  other  minor  has  the  same  order.  Hence,  (252)  can  be  written,  -  in  a  <  x  <  t  -,  as 
follows: 


(x,t, p) 


T^-n+l 

(-l)"-^^1  I,  .(VI 


-p  f  *{5>dc 


<f  g4’,(a)  n(a) 

^"“’(t)  n(t) 


v(x) 


(265 


where 


(u,(x)l . (u  ,(x)] 

I  ll"  1 

‘VV . lWlJ 


v(x) 


(L^tu,  )1 

‘WV 

•  •  • 

(L  (u.)J 

II  1 


Wi11 


[tr+2<Vin 


(L 


n(un-in 


(266) 


is  the  cofactor  of  the  minor  ±2  L^V)  I^^tv)  in  (252).  (265),  (258)  and  (254)  show 

that  for  t  <  x. 


G(x,t,p) 


(-1) 


vrn+1 


-p  /*  <^( oae 


±  ^(a)  n(a)  v(x) 

v>n_,(t)  n(t)  a 


(267) 


and  therefore 


/*  G(x,t,p)p  f ( t)dt 


T_.-n+2 

(-1)  r+1  *  r+1 


(q)  n(o)  ^ ~  f  *(e)dC  P(t,p) 
A 


(268) 


with 


F(t,p) 


tit) 

[*n_1(t)  n(t) J 


The  integral  in  the  right  hand  member  of  (268)  tends  to  zero  as  p  since  its 

integrand  tends  to  zero  in  the  interior  of  the  interval  and  remains  bounded  at  the 

Tr+l"n+2 

endpoints •  6  remains  bounded,  because  we  have  assumed  that 

Tr+1  <n'2  ‘ 

Consequently,  the  left  member  of  (268)  tends  to  zero,  as  p  ♦ 

b)  In  the  second  integral  of  (262)  we  have  t  >  x.  Hence,  in  the  determinant  in 
(252)  the  last  term  of  the  first  row  is  -V(x,  p).  Expanding  this  determinant  in  terms  of 
the  minors  formed  by  its  first  and  last  columns  we  see  that  this  time  the  minor  of  highest 
order  is  formed  by  the  two  first  rows  of  the  determinant.  The  value  of  this  minor  is 

±2  V ( x ,  p)  L^v)  . 

A  calculation  analogous  to  that  used  in  part  a)  of  this  section  leads  to 

JB  G  (x,t,  p)  p  f(t)dt  -  p""*2  [A'**  2]  (269) 

1 

where  the  constant  Q  is  the  limit  of  the  cofactor  of  the  minor  above.  Since  t  >  x  in 

the  integral  in  the  right  member  of  (269),  the  integrand  of  that  integral  tends  to  zero  in 

the  interior  of  the  interval  of  integration,  x  <  t  <  B,  as  p  approaches  infinity. 
Furthermore,  the  integrand  is  bounded  at  the  endpoints  of  the  interval  of  integration. 
Consequently,  the  integral  tends  to  zero.  On  the  other  hand,  our  assumption  m  >  0 

implies  n-2  >  0,  and  the  power  of  p  in  (269)  is  therefore  not  positive.  Hence,  the 

left  member  of  (269)  tends  to  zero,  as  p  ♦  «. 

This  completes  the  proof  of  (261), 
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47.  Ne  now  turn  to  the  asymptotic  calculation  of  i2(x,p),  (aaa  (260)).  No  nota  firat, 
on  expanding  tha  determinant  in  (253)  with  reapect  to  ita  firat  column,  that  H2(x,t,p) 
can  be  written,  in  a  <  x  <  0,  in  the  form 


where 


1,(7) 


Vx't'p)  *  7*0  [H<x*t)1 


(270) 


H(x,t) 


u,  (  X  ) 


W 


un_1(x)  g(x,t) 


l2(Vl>  V«> 


VV 


Ln(«  ,)  L(g) 
n  n- »  n 


(271) 


ie  the  limit  of  the  cofactor  of  the  element  L.j(V)  in  (253).  Substituting  (270)  in  (259) 
it  follows  that 


Gj(x,t,  p) 


(_t  )n  i  **<*>0 
P  <^( t >  *A 


and  this  formula,  together  with  (260)  givea 

v*,p)  *  <-un  fe  ti^if(t)dt  . 

*t)7i 


In  this  expression  we  may  pass  to  the  limit  under  the  integral  sign,  since  the  asymptotic 
expression  for  the  integrand  is  valid  in  the  whole  interval  a  <  t  <  0.  Thus  we  obtain, 
replacing  at  the  same  time  *t)  by  its  value  -bg(0. 


11m  Z2<x,p) 

p*w 


H(x,t)  f(t) 
A 


bQ(t) 


dt 


(272) 
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Now  we  combime  (272)  with  (261)  and  (260),  end  conclude  that  the  right  member  of  (272)  Is 
the  value  of  lim  z(x,t,p).  But  the  definition  of  H(x,t)  and  A  in  (271)  and  (255), 

p-MO 

respectively,  show  us  that 


n-1  H(x,t) 


(-1) 


is  exactly  the  Green's  Function  belonging  to  the  differential  expression 

1 


bQ(x) 


K(y) 


and  to  the  boundary  conditions  ^  ( i  ”  2,3,...,n). 

Since  these  boundary  conditions  do  not  involve  derivatives  of  higher  than  (n-2)nd 
order,  this  proves  that  z(x,p)  tends  to  a  solution  of  the  limiting  differential  equation, 
satisfying  all  boundary  conditions,  except  the  first  one. 


3 
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48.  We  summarize  the  results  of  thia  (  In  the  following  theorem 


Theorem  5i 

We  consider  the  problem  (N )  defined  by  the  differential  equation 

-  N(y)  +  M(y)  -  f(x)  232) 

P 

and  the  boundary  conditions  (108).  N(y)  and  M(y)  are  differential 
expressions  of  the  form  (103)  and  (104)  with 

n-m-1  . 

We  make  the  following  assumptions! 

(a)  Conditions  1*  -  5*,  7*  and  8*  of  the  Main  Theorem  (section  4) 
are  satisfied. 

(b)  b„(x)  <0  in  a  <  x  <  8. 

<e>  tr+f  <  n-1 . 

(d)  0  <  r  <  n. 

(e)  f(x)  is  lntegrable  in  a  <  x  <  0. 

Then  the  solution  z(x,p)  of  the  problem  (N)  tends  -  as  p  approaches 
infinity  -  to  a  solution  of  the  differential  equation 

M(y)  -  f(x) 

satisfying  all  the  boundary  conditions  (108)  except,  in  general,  the  first 
one,  L,,<y). 

A  strictly  analogous  theorem  holds  for  (x)  >  0. 
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Chapter  III 


SOME  RELATED  PROBLEMS 

§1.  An  Example  for  Boundary  Laver  Problems  In  Systems  of 
Differential  Equating. 

49.  A  great  number  of  unsolved  boundary  layer  problems  with  Important  applications  can  be 
formulated  for  systems  of  ordinary  differential  equations.  He  are  going  to  discuss  In  this 
section  a  very  elementary  example  In  order  to  give  an  Idea  of  the  boundary  layer  phenomena 
that  can  arise  for  systems! 

He  shall  discuss  the  system 


-1  „ 

p  u" 

v" 


“  a  u  t  b  v 
•  c  u  +  d  v 


(301) 


with  constant  a,  b,  c,  A,  assuming  that 

a  !<  0  .  (302) 

As  boundary  conditions  we  prescribe 

u( a)  ■  u  ,  v(  a)  -  v  ,  u(  B)  ■  u_,  v(B)  -  v„  (303) 

a  a  p  p 

where  u  ,  v  ,  u„,  v„  are  constants, 
a  a  8  8 

The  "limiting  problem",  obtained  by  setting  p  1  -  0  in  (301)  is  equivalent  to  the 
differential  equation  of  second  order 

v"  “  a  v  -  0  (304) 


where 


and  to  the  relation 


1305) 


(306) 


The  following  questions  ariset 

(a)  Do  the  solutions  U(x,p),  V(x,p)  of  (301)  and  (303)  converge  as  p  +  «? 

(b)  What  boundary  conditions  do  the  limit  functions  U(x)  and  V(x)  satisfy  if  they 
exist,  and  are  they  solutions  of  the  limiting  differential  equation  (304)? 

The  functions  U(x)  and  V(x)  cannot  be  expected  to  satisfy  all  four  boundary 
conditions  (303)  and  also  the  condition  (306),  for  the  prescribed  boundary  values  (303) 
will,  in  general,  not  satisfy  ths  condition  (306).  The  answer  to  these  questions  is 
supplied  by  the  following  Theorem! 


Theorem  6» 

Let  u  -  U(x,p),  v  -  V(x,  p)  be  solutions  of  the  system  of 
differential  equations 


p"V  -  a  u  ♦  b  v  >> 

v"»cu  +  dv-' 


(a,  b,  c,  d  constants)  satisfying  the  boundary  conditions 
u( a)  »  u#,  v( a)  -  va,  u(  B)  “  ug»  v< 


(301) 


(303) 


where  u  ,  u_,  v  ,  v_  are  constants.  Let  us  further  assumes 
a  B  a  6 

Assumption  1 :  a  f  0  . 

Assumption  2s  a,  b,  c,  d,  a  and  B  are  given  in  such  a  way 
that  the  differential  equation 


y" 


(307) 


with 

5  -  I*  51  (305) 

and  the  boundary  conditions 

y(  u)  »  0  ,  y<  B)  “  0 

can  be  satisfied  by  the  function  y(x)  =  0  only. 
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Then  we  can  make  the  following  statements: 

(A)  The  function  V(x,p)  converges,  as  p  +  ■»,  to  a 
function  V(x),  satisfying  the  differential  equation  (307) 
and  having  the  boundary  values 

V(  a)  -  v  ,  V(8>  -  v,  (308) 

a  8 

(B)  If  a  >  0,  then  the  function  U(x,p)  converges  in 
a  <  x  <  8,  as  p  +  to  the  function 

0(x)  -  -  |  V(x)  (309) 


(which,  of  course,  satisfies  the  differential  equation  (307)). 

(C)  If  a  <  0,  then  the  function  U(x,p)  does  not  converge 
but  remains  bounded,  as  p  ♦  ",  except  when  the  prescribed  boundary 
values  satisfy  the  condition 


a  u  +  b  v  «  0 
a  a 

(310) 

«  uB+  b  v8-  0  , 


In  which  case  statement  (B)  remains  true  for  a  <  0  also. 


50.  To  prove  this  theorem  we  start  from  the  observation  that  U(x,p)  and  V(x,p)  are 
both  solutions  of  the  differential  equation 

p-,{y<4*  -  <J  y"}  -ay'tiy-O  .  (311) 

In  fact,  if  we  multiply  the  second  differential  equation  of  (301)  by  b  and  substitute 
into  it  the  expressions 

^  y  p  P  ^  U.  —  ^  |  J  ^  y,  ^  ^  ^  p  (  ^  ^  ^  p  W 

obtained  from  the  first  differential  equation,  we  find  the  following  differential  equation 


in  u  alone 


I 


u  -<pa  +  d)u"tpSu-0  , 

which  is  equivalent  with  (311).  Similarly,  it  can  be  shown  that  (311)  is  satisfied  by 
V(x,p),  by  eliminating  u  from  the  differential  equations  (301).  It  is  also  easily  seen 
that  the  boundary  conditions  satisfied  by  U(x,p)  and  V(x,p),  considered  as  solutions  of 
(311),  are,  respectively. 


for  0(x,p) 


for  V(x,p) 


y(a)  «  u 


y(B)  -  u„ 


y“(a)  •  (a  u  ♦  b  v  )  p 
a  a 


y"(0)  -  (a  ufl  ♦  b  vfl)p 


(312a) 


y(a)  -v 


y(B)  “  v. 


y"  ( a)  -  c  u  +  d  v 
a  a 


y-($)  -  c  U(J  ♦  d  vg  . 


(312b) 


The  boundary  conditions  for  V(x,p)  are  of  the  type  considered  in  our  Main  Theorem 

(chapter  I,  section  3).  Applying  the  Main  Theorem  for  these  boundry  conditions  and  for 

n-m  -  2,  bg ( x)  -  -a,  we  see  that  for  a  >  0  as  well  as  for  a  <  0  the  function  V(x,p) 

converges  to  that  solution  of  the  limiting  differential  equation  (307)  which  satisfies  the 

boundary  conditions  y (a)  “  v  ,  y ( 6)  -  v  .  This  completes  the  proof  for  statement  (A). 

a  B 


51.  For  U(x,p)  a  special  calculation  is  necessary,  since  the  boundary  conditions  for 
U(x,p)  contain  p,  a  case  not  considered  in  our  Main  Theorem.  Our  principal  tool  in 
chapter  I,  the  asymptotic  representation  of  a  fundamental  system  of  solutions  of  the  given 
differential  equation,  can  be  applied  to  the  differential  equation  (311)  and  yields  then 


(see  theorem  1,  section  6) 


U1  (x, p)  -  e 


oVa(x-a) 


-o^a(  x-a) 


Uj(x,p)  -  e 


Uj(x.p)  -  (A(x)l 


a  <x,p)  •  (B<  x)) 


Here  the  functions  A( x) ,  B(x)  form  a  fundamental  system  of  solutions  of  the  limiting 
differential  equation  (307)  and 


a  m  /p  • 

With  these  functions  we  repeat  the  general  reasoning  of  chapter  X  for  this  special  case. 

We  first  re-write  the  boundary  conditions  (312a)  in  a  form  somewhat  more  convenient  for  our 
purpose : 


-1 

P 


y"(  6) 


1.  “au„+bv„ 
1  8  8 


-1 

P 


y"(  a) 


-  a  u  +  b  v 
3  a  a 


y(  8)  -  t2  -  ug 


y(  a)  *  t.  »  u 

4  a 


Let  us  further  introduce  the  abbreviations 

e^(  0-a)  ,  ^ 


.  e'^(  a'a)  -  W2  . 


(314) 


(315) 


Plrst  case:  a  >  0. 

Defining  A(  p)  as  in  (134)  we  find 


(a)Wt 

faJWj 

f'tA’lB)) 

A(  p)  - 

(1)W1 

(1!W2 

[A(  8)J 

[8(  8)) 

(316) 

[a] 

(a) 

p_1tA-(o)] 

p'’[B"(a)] 

(11 

m 

[A  (  a)  1 

[B  (  a)  ] 

• 

We  expand  this  determinant  in 

terms  of 

its  last  two  rows: 

A(p>  -  -[a2] 

At  8) 

Ato) 

B(  B) 

w  » 

1 

B(  a) 

-[a2]D  Wt 

(317) 

with 


D 


A(6)  B(8> 

A(  a)  B(  a) 


In  consequence  of  assumption  2  of  the  theorem  to  be  proved  we  have 


(318) 
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I 


Similarly,  vie  find 


41 

[a]W2 

p_1[v<e>i 

p"1[B*(0)l 

*2 

(i]w2 

[a<  en 

[B(  6)1 

*3 

[a] 

p_1[A-(a)l 

p'1  [B*(a)l 

*4 

HI 

[A(  ail 

[B(  a)  1 

taiw. 

S 

p"1  [As  (0)1 

P_1  [B«(0)1 

inn, 

*2 

[A(  0)1 

[B(  0)] 

[a] 

*3 

p  1 [A* ( a)  1 

P_1 [B* ( a) ) 

[11 

*4 

(A<  a)) 

[B(  a)l 

-ta  11  D] 


-  -[a  *3  D] 


In  the  determinant 


[alW, 

[aJW2 

*1 

p"1[B"(0)l 

[ilw1 

li]i»2 

*2 

[B(  0)] 

[a] 

[a] 

*3 

p”1 [B* ( a)  ] 

[11 

Ml 

*4 

[B(  a)  1 

we  subtract  the  a”  fold  of  the  first  roe  from  the  second  row  and  the  a  fold  of  the 
third  row  from  the  fourth,  thus  obtaining 


[alW1 

[a)W2 

*1 

P_1[B*(0)1 

[01W1 

[0]w2 

b 

-  —  V 

a  a 

[B(  0)  1 

[a] 

[al 

S 

p_1[B"(a)l 

[0] 

[0] 

b 

-  —  V 

a  a 

[3(0)1 

Here  we  have  used  the  fact  that,  in  consequence  of  013) 


Expanding  (321)  we  find 


Similarly,  we  prove  that 


Aj(  0) 


b 

a 


B<  S) 

B(  a) 


(a2)W, 


A<8> 

A(  a) 


[a2]W. 


(322) 


(323) 


(324) 


52.  From  (313),  (317),  (320)  and  (321)  it  follows  that 


lim 

p-M» 


A,(p) 

A(p) 


O^X.p) 


0 


lim 

p-*°o 


*2(P) 

A(p) 


«2<x,p) 


0 


in  o  <  x  <  0  . 


(323)  and  (324)  show  then  that 

U(x)  -  lim  0(x,p) 

P*» 

Is  that  solution  of  (307)  which  assumes  the  boundary  values 

U(c)  -  -  ~  v  ,  5(8)  -  7  v  . 

a  a  a  b 

This  proves  statement  (B)  of  theorem  6. 

Second  case;  a  <  0  . 

If  a  <  0,  then  the  absolute  values  of  w,  and  W2  oscillate  with  increasing  p 
without  tending  to  a  limit.  In  this  case  we  obtain  from  (316)  the  asympotlc  expression 
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Co 


(325) 


A2)  (W,  -  V»2) 

instead  of  <317) .  For  A^p)  and  AjJp)  we  have 

A, ( p)  -  -CO  a]  <Jt1  -  <3  <*2) 

Aj(p)  -  -CO  a]  (t3  »1  -  t3) 

and  for  A3<p)  and  A^(p) 


Aj  ( p) 


-C«21 


b 

“  a  VB 


b 

—  v 
a  a 


B(  6) 

B(  a) 


(W,  -  w2) 


a4(p)  -  -[a  1 


A(  0) 


A(a) 


b 

a  VB 


b 

a  vo 


<W. 


V 


53.  From  these  expressions  It  follows  immediately  that 


4  A  (p) 

U(x,p)  -  —  0l(x,p) 


does  not  converge  in  this  case.  For 

77^7  <&3<P>  ♦  &4(P)  V*'p)  > 

converges  to  the  same  solution  of  (307)  as  in  the  case  a  <  0,  while  the 


J,  /  •  u  um.  i/T*l(  x-a)  VVS  -io/TIl( 

A  <41U1  +  42°2  a(W,^2)  *  +  *^V^T) 


(326) 

(327) 

(328) 

(329) 


expression 

x-a) 
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does  not  converge  unless  4^  “  t  ••  0,  in  which  case  it  is  identically  zero. 

This  completes  the  proof  of  theorem  6. 

If  assumption  2  of  theorem  6  is  not  satisfied,  then  our  reasoning  does  not  hold  any 
more.  In  that  case  it  would  be  necessary  to  take  into  consideration  also  the  second  terms 
of  the  asymptotic  series  used,  in  order  to  find  the  order  of  magnitude  of  4(p). 

Assumption  2  is  easily  seen  to  be  equivalent,  in  this  case,  with 

Si>0 


and 


J-  (f5-a)  +  Nrt  when  ~  <  0 


where  N  is  any  positive  or  negative  integer. 

The  case  a  »  0  could  be  easily  treated  by  the  same  method. 

The  boundary  layers  in  this  §  occur  only  for  the  function  U(x,p),  which  with 
increasing  p  tends  to  a  function  which  does  not  have  the  prescribed  boundary  values, 
except,  when  these  boundary  values  satisfy  the  condition  (310). 

A  more  adequate  and  general  treatment  of  boundary  layer  problems  in  systems  of 
ordinary  linear  differential  equations  could  probably  be  based  on  the  asymptotic  solution 
of  linear  systems  as  developed  by  Langer  and  G.  D.  Birkhoff  [51.  The  assumptions  of  that 
theory  would,  however,  have  to  be  generalized  for  this  purpose. 


$2.  An  Example  for  Boundary  Layer  Problems  with 
Singularities  in  the  Interior. 

54.  Introduction!  If  the  assumption  6*  of  the  Main  Theorem  in  chapter  I  is  dropped,  i.e. 
if  we  admit  zeros  of  bQ(x)  in  a  <  x  <  B,  then  our  whole  theory  becomes  Invalid.  For 
the  zeros  of  bg(x)  are  usually  singularities  of  the  limiting  differential  equation 
M(y)  -  0,  and  the  theory  of  the  asymptotic  solution  of  differential  equations,  which  was 
our  main  tool,  fails  in  this  case. 
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The  general  treatment  of  boundary  layer  problems  In  thla  caee  would  probably  require 
an  entirely  new  approach.  But  It  la  already  Interesting  to  Investigate  a  very  simple 
special  problem  of  this  type  in  which  the  differential  equation  can  be  solved  explicitly. 
We  shall  discuss  the  boundary  layer  problem  of  the  differential  equation 

-  y"  ♦  b(x)y'  -  f(x)  (330) 

P 

with  the  boundary  conditions 

y(o)  -  i2  ,  y(S)  «  i,  .  (331) 

We  make  the  following  assumptions : 

1.  b(x)  is  regular  analytic  in  a  <  x  <  0  . 

2.  f(x)  is  regular  analytic  in  a  <  x  <  6  . 

3.  b(x)  has  a  zero  at  the  interior  point  x  «  r  of  the  interval  a  <  x  <  B.  But 

b' (r)  0.  and  b(x)  does  not  have  any  other  roots  in  a  <  x  <  8*  (This  number  r  has, 

of  course,  nothing  to  do  with  the  number  r  used  in  the  first  two  chapters.) 

Assumptions  1  and  2  are  by  no  means  essential.  We  introduce  them  only  in  order  to 
simplify  our  reasoning. 


55.  The  solution  of  the  boundary  value  problem  for  the  differential  equation  (330):  To 
simplify  the  calculations  it  shall  be  assumed  that  the  boundary  values  are 

y(  a)  -  y(  6)  -  0  . 

The  case  of  non-hosK>geneous  boundary  conditions,  which  does  not  add  any  new  features  to  the 
problem,  la  discussed  in  section  59. 

The  general  solution  of  (330),  as  obtained  by  elementary  methods  can  be  written  in  the 

form 


(J(x,p)  -  J*  d?  /*  p  f(n)e'p(M5)-Mn))dn  +  c1  J*  e_pMe>de  +  c2 


where 


A(x)  -  b(x)  dx  (333) 

and  X,  y,  v,  ,  c ara  arbitrary  constants,  which  ara  rastrlctad  by  tha  prescribed 
boundary  conditions.  Tha  five  constants  ara,  of  coursa,  not  assantial.  In  reality  (332) 
depends  only  on  two  essential  parameters,  so  that  thraa  of  the  five  constants  can  be  chosen 
arbitrarily.  In  order  to  obtain  a  form  of  tha  solution  suitable  for  the  boundary  layer 
problem.  It  Is  convenient  to  set  u  “  1,  while  the  choice  of  X  and  v  shall  be  left 
undecided  for  the  moment.  For  typographical  reasons  It  Is  useful  to  Introduce  the 
abbreviations 


ps  “  £  de  ^S  0  f<n)e"pU<C)"x(n))dn  (334) 

Q*  -  £  e‘pX(C)dC  .  (335) 

Then  (332)  can  be  written 

0<x,p)  -  P*  ♦  CjQ*  ♦  c2  ,  (336) 

and  substitution  of  tha  values  x  ”  a  and  x  »  8  into  (336)  leads  to  the  two  linear 
algebraic  equations 


0  *  p  ®  +  c  0  ®  c 
X  1  X  2 

0  ■  +  e  +  e 

X  1WX  2 


for  c1  and  c2.  Calculating  c2  from  these  equations  and  using  the  fact  that 

Q®  -  q“  »  Q®,  one  finds  that 
A  A  a 


pxB®x  -  p?> 


(337) 
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investigation  can  ba  obtained  immediately  fro*  (337)  by  tha  following  conaidarationai 
Substitution  of  X  for  x  in  (336)  shows  that  c2  -  u(X).  As  X  was  arbitrary,  this  ia 
true  for  any  value  of  X,  so  that  (337)  can  ba  regarded  aa  the  desired  solution  of  (330) 
with  X  instead  of  x  as  independent  variable.  Writing  x  for  X,  the  solution  of  the 
boundary  value  probla*  is  therefore  obtained  in  the  fora 


,  °*„a  «*„. 

a  a 


56.  The  asyaptotic  value  of  /  F(  x) e  dx  for  large  pi 

The  solution  of  (336)  is  composed  of  integrals  of  the  for* 

P(x)epl*x>dx  .  (339 

It  la  therefore  laportant  to  have  asymptotic  expressions  for  such  integrals  for  large 
values  of  p.  The  following  theorem,  a  proof  of  which  can  e.g.  be  found  in  a  paper  by 
O.  Perron  [8],  will  be  the  chief  tool  of  the  subsequent  investigations. 

Theory* i  If  F(x)  and  <4x)  are  regular  analytic  in  s  <  x  <  t,  if 


C  «  0,  for  x  •  R,  a  <  R  <  t 
*<x)  < 

v  <  0,  for  x  R,  s  <  x  <  t 


<p(x)  »  (x-R)p(gQ  +  g^(x-R)  +  •••) 


is  the  Taylor  series  of  (x)  around  x«R,  then 


ft  F(x)epl*x)dx  -  if  r(R)r4>  - W)  p_Vp  , 

JS  P  P  .  |VP 


P  l90l 


where  the  brackets  "(  ]*  have  the  meaning  defined  in  section  5.  (The  number  p  here  has 
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where  tha  brackets  "t  1"  hava  tha  meaning  defined  in  section  5.  (The  number  p  hare  hae 
nothing  to  do  with  the  number  p  of  chapter  I.)  If  R  "  s  or  R-t,  tha  same  is  true, 
but  with  the  factor  2  in  (340)  missing. 

This  theorem  can  easily  be  generalized  so  as  to  include  also  the  case  ^(x)  >  0: 
Theorem i  If  P(x)  and  *(x)  are  regular  analytic  in  s  <  x  <  t,  if  *(R)  is  the 
maximum  of  *>(x)  in  s  <  x  <  t,  where  s  <  R  <  t,  and  if 

*>(x)  -  *<R)  +  (x-R)p(gQ  +  (x-R)g1  *  ...) 

is  the  Taylor  series  of  ^ (x)  around  x  “  R,  then 

f*  F(x).^lx)dx  -  [i  P(R)r<i)  — L— 1  p-1/P  .*<*>  .  (341) 

p  p  lg0l1/p 

If  R  »  s  or  R  -  t,  the  same  is  true  with  the  factor  2  in  (341)  missing. 

Proof i  The  integral 

£  P(x)eP<*<X)-*(R))dx 

satisfies  the  assumptions  of  Perron's  theorem  with  (x)  -  v  (R)  instead  of  i p(x).  Since 

*(x)  -  *(R)  -  (x-R)p(g0  +  (x-R)g1  ♦  ...)  , 

(341)  follows  immediately,  if  (340),  applied  for  the  exponent  <Mx)  -  f  (R)  is  multiplied 
on  both  sides  by  ep  . 

57.  Passage  to  the  limit  in  (338),  if 


<  0 

for 

a 

<  x  <  r 

b(x)  -  0 

for 

X 

-  r 

>  0 

for 

r 

<  x  <  $  • 
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The  asymptotic  calculation  of  the  integrals  in  (336)  is  easiest  when  the  constant  v 
in  the  definition  of  A(x),  formula  (333),  is  chosen  equal  to  r.  If  thia  is  done,  the 
function  A(x)  satisfies  the  following  conditions  > 


A(x)  >  0,  x  j*  r 
A(x)  «  0,  x  -  r 

A( x)  is  monotonic  increasing  for  x  >  r 
A(x)  is  monotonic  decreasing  for  x  <  r 


(342) 


M  x)  has  therefore  the  shape  indicated  by  the  figure  below. 


The  passage  to  the  limit,  as  p  *  «•,  in  (336)  leads  to  different  results,  according  as 
x  is  less,  equal  or  greater  than  r. 

Case  a) .  x  >  r. 

x  B 

Application  of  formula  (341)  to  Q  and  Q  show  that  for  x  >  r  both  integrals 

B  a 

have  the  same  asymptotic  value,  hence 


lim  — ^  »  -  lim  — 2  -  -1  •  (343) 

P+*  8„  P**° 

a  a 


a 

To  the  inner  Integral  of  formula  (340)  can  be  applied.  For,  in  this  case, 

r  <  x  <  n  <  5  and,  in  this  range,  —  < A< O  -  A(  n) )  as  function  of  n  assumes  its  maximum 
value  0  for  n  •  £,  on  account  of  (342).  As 


|n(-A(  5)  +  A(  n) 


-b(  n) 


is,  by  assumption,  not  zero,  the  number  p  of  (340)  is  here  equal  to  1,  hence 
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J*  p  f(n>e“ptMO"M’'Mdii-  (f(£) 


and  therefore 


lim  P® 

P-Mi  X 


f®  iiii 

jx  b<£) 


*5 


(344) 


The  expressions  in  brackets  in  (340)  and  (341)  can  be  sure  to  be  different  from  zero, 
if  F(x)  r1  0  in  a  <  x  <  6*  In  order  to  avoid  too  lengthy  formulas,  this  additional 
assumption  shall  temporarily  be  made.  It  is,  however,  by  no  means  essential,  and  it  will 
be  shown  later  how  to  proceed  without  it. 

The  letter  E  in  the  subsequent  formulas  shall  be  used  to  indicate  non-vanishing 
constants.  Note  that  the  same  letter  E  will  be  used  for  different  constants. 

One  finds  immediately,  by  application  of  (341), 


Q 


B 

x 


-pMx) 


(345) 


in 

Q8-  /».-**  *>dC 
a  1  a 

the  exponent  reaches  its  maximum  0  for  £  “  r.  As,  by  assumption,  b(r)  *  0,  but 
b'(r)  f*  0,  p  is  equal  to  2  in  this  case  and 

Q?  -  »I  7T  •  (346) 

a  v  p 

In  order  to  find  the  order  of  magnitude  of  p“,  consider  that  a  <  £  <  n  <  x  and  r  <  x 
in  the  exponent  -(A(£)  -  A(r\))  occurring  in  P°.  Hence,  the  maximum  of  -(A(  £)  -  A(  n) ) 
as  function  of  n  for  fixed  £  is 


-<A(£) 

-  A(  £))  -  0, 

for 

A(£> 

>  A(x) 

”( A(  £) 

-  A(x) )  >  0, 

for 

A(  £) 

<  A(x) 

The  asymptotic  value  of  the  inner  integral  of  P*  is 
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/epf(n).-p(A(5,"MT’))4n 


IB] 


,  for  A(  C)  >  A(x) 


[B}*-ptM5>"Ml0)  ,  for  A(£)  <  A(x) 


The  contribution  to  of  that  part  of  the  interval  a  <  £  <  x  for  which  A(  £)  >  A(x), 

(if  it  exists) ,  can  be  neglected  in  comparison  with  the  part  where  A( £)  <  A(x).  As 


/“  .-P(A(5)-A(x))^  .  m  ^ 


1  pA(x) 


one  haa  therefore 


(B1 


1  pA(  X  ) 


75“ 


(347) 


From  (345),  (346)  and  (347)  It  follow*,  finally,  that 


P“  -  IB]  l 

Q6  x  P 

a 


and  therefore,  using  (344), 


Um  W(x, p)  -  -  /®  dC,  for  x  >  r 

P-H* 


(348) 


Caae  b) .  x  <  r . 

A  consideration  analogous  to  that  used  for  case  a)  leads  to 


11m  U(x,p)  -  /*  £7^  d£,  for  x  <  r  . 

”(5) 


(349) 


See  also  case  c)  and  d) ,  section  58. 

Case  c) .  x  ■  r. 

Is  regular  analytic  at  x  «  r,  then  one  shows  Immediately,  by  a  similar 
that 


if  £rr 

b(x) 

consideration. 


2  b(f}  ^  *  1 6  bff)  d5*'  for  x”r 
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(The  factor  j  is  due  to  the  Q's  in  (338).)  in  other  words,  U(x, p)  tends  in  this  case 
to  the  arithmetic  mean  of  the  two  limits  at  x  -  r. 


58.  Passage  to  the  limit  in  (338),  if 

C  >  0,  for  a  <  x  <  r 

b(x)  \  -  0,  for  x  »  r 

\  <  0,  for  r  <  x  <  6  . 

In  order  to  operate  as  much  as  possible  with  positive  quantities  it  is  convenient  to 
set  now 


A(x)  -  J*  b(x)  dx  . 


(350) 


Then  A(x)  satisfies  the  conditions  (342).  If  in  the  definition  of  Pfc  and  o'"  the  sign 

8  s 

of  the  exponents  is  changed,  U(x,p)  can  again  be  written  in  the  form  (338). 

In  addition  to  the  distinction  between  the  cases  x  >  r  and  x  <  r,  the  relative 

size  of  A(  a)  and  A( 8)  plays  now  a  part  in  the  proof. 

Case  a],  x  >  r,  A(8)  >  A(a)  . 

I*t  s  >  r  be  the  value  for  x  for  which 


A(s)  -  A(  a) 


Then 


(E]  -  epA(X)  ,  x  >  s 


[E]  ’  e^3’  ,  x  < 

0 


(351) 


and 


Q8  -  [E]  -  epA<6)  . 
a  0 

o 

The  relative  positions  of  x,  £  end  n  in  P®  are  indicated  in  the  figure  below. 


(352) 
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One  sees  that  max(A(£)  -  A(r>))  -  A<?)  -  A(x)  for  fixed  i.  Therefore 


/J  p  f(n>.pUU)-Mn))dn-  [■]  eptMe)-A(x)) 

rS  bP<A<5)-A<x))  2  eP(A(  S)-A(X)) 

■*x  P 


and 

_8  1  p(  A(  8)-A(x) ) 

Px  -  [E]  -  e 

Furthermore 

For  pa  consider  again  the  relative  position  of  x,  £  and  n*  If  t  > 
Integral  in  P®  remains  finite,  as  p  ♦  »,  hence  only  the  case  £  <  r  has 
considered#  In  that  case 

max(A{£)  -  A(n))  -  A(C)  -  Mr)  -  M  C) 

and 


r  p  f(n) 


p(a<  e)-R(  h>)d 


tE]  /p  e 


pM£) 


As 


one  has 


f°  epA<t)dt  - 
Jx 


rtE]  -;epA(x) 

\[I]  -1pepA(a) 


x  >  s 


x  <  3 


(353) 

(354) 

the  inner 
to  be 
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V. 


1  pA(  x) 

[El  e  ,  x  >  s 


[E]  e  ,  x  <  a 


Substituting  formulas  (351)  -  (338)  ona  obtains 


U(x,p) 


tE1  i 


[El  A.  apA(X) 

*P 


(El  -  Bp<A<a)-M*)>  t_,  1  pA(a) 

1  J  p  "  lE)  /p  8 


lim  U(x,p)  -  ±»  . 

p-ooo 


Case  b).  x  >  r,  A(B)  <  >(a)  . 

The  reasoning  is  of  the  same  type  as  in  case  a) ,  only  the  orders 
terms  change.  Let  s  <  r  be  the  value  of  x  for  which 


One  finds 


*(a)  -  A(  8) 


—f,  -  [El  . 

Q8 

a 


Por  px  the  asyraPtotlc  formula  (353)  holds  unchanged.  Furthermore 


QS  -  [El  i  epA(B) 

X  p 


Q8  -  [El  1  epA(a) 

a  p 


The  relative  positions  of  x,  £  and  n  in  the  inner  integral  of 

x 

figure  below.  Only  the  case  £  <  r  has  to  be  considered,  as  for  £  : 


(355) 

x  >  s 

X  <  s 

(356) 

of  magnitude  of  the 

(357) 

(358) 

(359) 

can  be  seen  in  the 
r  the  integral 


tends  to  zero. 
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JJ 


1 


For  e  <  r 


jj  P*(rt  .<*M»-w»d,.  m  /p.^n 


/a.(A(5)d5-  m  ip <*<“> 
x  p 


The  expression*  (353)  and  (357)  -  (360),  Inserted  In  (338)  lead  to 


U(*.«)  -  -  m  Ae.^W 

O  *0 


Therefore 


11a  U(x,p)  -  *••  . 

P*~ 


Ca*e  c)  and  d).  x  <  r. 


One  might  repeat  the  preceding  argument*  in  analogous  fora  and  with  the  same  final 
result  for  x  <  r.  Instead,  one  can  also  proceed  a*  follows i  The  transformation 


t  m  r  -  X 


changes  the  differential  equation  (330)  into 


+  *(i)  £■ ?(z) 

02 


b(*)  »  -b(r-s) 


f(*)  •  f(r-s) 


-92 


The  boundary  conditions  yta)  *  y(  B)  “  0  are  transformed  into 


As 


y(r-  B)  “  y(  r-a)  -  0  . 


f  * 

0, 

for 

r-  B  <  * 

<  0 

b(z)  \  » 

0, 

for 

z 

-  0 

l  < 

0, 

for 

0  <  z  < 

r-a 

The  results  of  the  preceding  section  can  be  applied  to  (362)  and  lead  to 

11m  U(z,p)  •  ±«*  ,  for  0  <  z  <  r-a  • 

P-m* 

But  this  is  equivalent  with 

lim  U ( x,  p)  -  ±“  ,  for  a  <  x  <  r  . 

P-h. 

This  method  might  also  have  been  used  to  prove  formula  (349). 


59.  Some  minor  generalizations: 

1)  If  f(x)  is  allowed  to  have  roots  in  a  <  x  <  B»  the  asymptotic  values  used  for 

p  and  Pp  are  not  correct  for  these  values  of  x  which  are  roots  of  f(x).  If,  in 
x  x 

particular,  f(r)  ”  0,  some  of  the  expressions  would  even  be  incorrect  for  all  x.  The 
reason  is  that  f(x)  or  f(r)  appear  in  the  constant  factors  occurring  in  these 
expressions.  But  in  that  case  the  calculation  could  be  carried  through  with  a  slightly 
more  general  form  of  formula  (340)  also  contained  in  Perron's  general  formula  in  the  paper 
(8]  quoted  above.  The  result  is  again  the  same. 

2)  Non-homogeneous  boundary  conditions:  If  to  (338)  is  added  a  solution  of  the 
homogeneous  differential  equation 

-  y"  ♦  b( x)y'  -  0  (363) 

P 

satisfying  the  non-homogeneous  boundary  conditions 

y(o)  -  tj,  y(S>  -  l,  (364) 

one  obtains  the  solution  of  (330)  satisfying  the  boundary  conditions  (364). 
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I 


4 

Mr 

* 


As  tha  transformation 


y  -  y  -  t. 


changes  (330)  into  a  differential  aquation  of  the  same  type  and  transforms  the  non- 
homogeneous  left  hand  boundary  condition  Into  a  homogeneous  one,  It  Is  no  loss  of 


generality  to  assume  that 


s-° 


The  general  solution  of  (363)  can  be  written  In  the  fora 


U(x,p)  •  c(  Q*  +  c2 


The  given  boundary  conditions  lead  then  to 


°1  “  li  '  c2-°  ' 

wa 


0(x,p)  -  —  t  . 

Q®  1 


To  (365)  the  methods  of  sections  57  and  58  can  be  immediately  applied  with  the  following 


result i 


r  <  0,  x  <  r 
l  >  0,  x  >  r  , 


0  for  x  <  r 


f  1.  for  x  >  r 
lim  tJ(x,p)  “  < 

p-M«  v.  0  for  x  <  r 

(>  Q,  for  x  <  r 
b).  b(x)  ( 

l  <  0,  for  x  >  r  , 


a.)  A(  B)  »  A(a)  , 


U(x,p)  - 


[«]  .P<*(x>-A(B)\  x>  , 
[«1  .p<*<«>-a(B))#  x  <  , 


. «  •  . 


M  0)  <  M  a) 


6.) 

lia  0(x,p)  -  t1  . 

P*** 


60.  Theorem  7i 

Given  the  differential  equation 

iy*  ♦  b(x)y'  -  f(x)  (330) 

P 

where  b(x)  and  f(x)  are  regular  analytic  functione  in  the 
interval  a  <  x  <  0  and  b(x)  has  exactly  one  root  x  -  r 
in  the  interior  of  the  interval,  while  b' (r)  4  0.  Then  the 
behavior  for  great  valuea  of  p  of  the  aolution  of  (330)  which 
aatisfies  the  boundary  conditiona 

y( a)  -  t2.  y< 0)  -  i,  (331) 

depends  essentially  on  the  shape  of  b(x): 

C  <  0  for  a  <  *  <  r 
1).  If  b(x)  < 

'  >  0  for  r  <  x  <  0  , 

then  the  solution  n(x,p)  converges  with  Increasing  p  in 
the  whole  Interval  a  <  x  <  0,  except  possibly  at  x  *  r,  the 
limiting  function  being  composed  of  the  two  solutions  of  the 
limiting  differential  equation  of  the  first  order 

b(x)y'  -  f(x) 

satisfying  one  of  the  two  prescribed  boundary  conditions.  If 
these  functions  are  bounded  at  x  «  r,  then  the  solution  of 
(330)  converges  at  x  »  r  to  the  arithmetic  mean  of  the  two 
limits  at  this  point. 

C  >  0  for  a  <  x  <  r 
2).  If  b(x)  < 

v  <  0  for  r  <  x  <  0  , 
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* 


then  D(*,o)  dlvargaa  with  lncraaalng  p  at  all  pointa  of  tha 
lntarval. 
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Appendix 

A  SHORT  REPORT  ON  THE  ASYMPTOTIC  SOLUTION  OP 
ORDINARY  DIFFERENTIAL  EQUATIONS  INVOLVING  A  PARAMETER 

61.  The  aaln  mathematical  tool  used  In  thia  paper  la  the  theory  of  the  asymptotic  solution 

of  differential  equatins  Involving  a  parameter  P  for  large  values  of  this  parameter.  The 

most  important  results  of  this  theory  are  contained  in  papers  by  G.  □.  Birkhoff  [ 1 ] , 

Noaillon  [2],  Tamarkin  13] ,  [4],  Langer  [5],  and  Trjitsinsky  [61.  The  asymptotic 

developments  used  by  Birkhoff  and  Tamarkin.  although  of  a  very  general  character,  do  not 

apply  to  the  particular  differential  equations  of  this  investigation,  because  they  assume 

that  a  certain  "characteristic*  algebraic  equation  formed  with  the  coefficients  of  the 

differential  equation  has  no  multiple  roots,  an  assumption  not  satisfied  in  our  case. 

It  would  probably  not  be  difficult  to  modify  the  methods  used  by  Birkhoff  and  Tamarkin 

in  such  a  way  that  they  cover  our  case.  But  this  is  not  necessary,  since  the  type  of 

differential  equations  considered  by  Noaillon  and  Trjitsinsky  includes  the  differential 
equation  (101). 

A  complete  proof  of  the  main  theorem  of  Noaillon  and  Trjitsinsky  would  be  beyond  the 
scope  of  this  investigation,  even  if  we  restricted  ourselves  to  the  special  case  in  which 
we  are  interested.  We  intend  here  only  to  give  a  summarised  report  on  the  methods  of  this 
theory  and  to  show  how,  assuming  the  theorems  proved,  the  asymptotic  expressions  of  theorem 
1  can  be  obtained  in  our  case. 

The  theory  consists  of  two  parts.  In  the  first  part,  which  may  be  called  the  "formal 
part",  the  "exact"  differential  equation 

L(y,p)  -  0  (401) 

which  is  essentially  equivalent  to  the  differential  equation 

— ^lp)-  -  0  (402) 

y 

is  replaced  by  the  "asymptotic"  differential  equation 

L(y,P)  „ 

Y  j 


0 


(403) 


Her*  the  symbol  "j"  has  the  following  meaning] 


'1 

•  * 

Jt 


f(x,p)  J  g(x,p)  , 

where  j  Is  a  real  number,  stands  for 

f(x,P)  -  g(x,P)  -  . 

p"3 

As  previously,  E(x,P)  Is  a  function  such  that  there  is  a  positive  real  number  R  so  that 
|R(x,P)|  is  uniformly  bounded  for  a  <  x  <  0  and  |p|  >  R.  Usually,  f(x,P)  will  be 
regarded  as  asymtotically  equal  to  g(x,P)  only  if  j  is  positive.  But  sometimes,  e.g. 
in  the  case  of  the  differential  equation  (101),  a  solution  of  (403)  with  a  negative  j 
will  be  an  asymtotic  approximation  in  the  ordinary  sense  of  the  exact  diferentlal  equation 
(401). 

It  is  then  shown  that  under  very  general  assumptions  a  function  Y(x,p)  can  be 
constructed  which  satisfies  the  condition  (403),  provided  the  number  j  is  not  too  large. 

If  certain  differentiability  conditions  are  satisfied,  j  may  have  an  arbitrarily 
large  value.  This  is  the  case  treated  in  detail  by  Noalllon  and  Trjitxinsky.  The  case  in 
which  there  is  an  upper  limit  for  j  is  only  mentioned  occasionally  by  these  authors.  But 
since  we  are  only  interested  in  the  first  terms  of  the  resulting  asymptotic  expansions,  it 
is  unnecessary  to  assume  indefinite  differentiability  of  the  coefficients  of  the 
differential  equation.  This  assumption  is  required  only,  if  we  are  interested  in  the 
unlimited  asymptotic  expansion.  Going  over  Noaillon's  proof  it  is  easily  seen  that 
assumption  4*  of  Noaillon's  theorem  in  section  62  below  is  sufficient  to  guarantee  the 
existence  of  the  first  term  of  the  asymptotic  solution. 

In  the  second,  the  "functional''  part  it  is  shown  that  the  solutions  of  the  asymptotic 
differential  equation  (403)  are  asymptotically  equal  to  the  solutions  of  the  exact 
differential  equation  (401). 

Essential  for  our  application  is  furthermore  the  result  that  the  derivatives  of  these 
asymptotic  solutions  of  (401)  are  asymptotically  equal  to  the  derivatives  of  the 
corresponding  exact  solutions. 
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Th«  complete  etateaent  of  the  result*  of  Bouillon's  paper,  as  far  as  they  are 


Important  for  our  purpose,  follows: 


62.  Soalllon's  Theorem. 

Part  Ii  Given  Is  the  differential  equation 

WfiP)  “  l  VAx,  p)y<n"1)  -  0  (404) 

1-0 

satisfying  the  following  conditions: 

1*.  p  is  a  real  positive  parameter. 

2*.  x  is  a  real  variable 

3«.  In  the  domain  a  <  x  <  8,  p  >  R  <  a,  8,  R  constants)  the  coefficients  P^fx, p)  can 
be  expanded  in  convergent  series  of  the  fora 

H  — 

PA(x,p)  -  p  1  l  B  (x)p“*  .  (405) 

s-0 

(The  are  positive  Integers.) 

4".  The  functions  B1#  have  at  least  n  continuous  derivatives  in  a  <  x  <  8. 

5«,  The  coefficients  Bao(x)  in  (405)  does  not  vanish  in  any  point  of  the  interval 
a  <  x  <  8* 

To  these  conditions  1*  -  5*  a  sixth  assumption  has  to  be  added,  which  can  be  most 
easily  defined  in  the  course  of  the  construction  of  the  asymptotic  solution. 

If  these  conditions  are  satisfied,  then  there  can  be  constructed  solutions  Y(x,p)  of 
(403),  each  of  which  can  be  Witten  in  the  form 

Y(x,p)  -  T*u  (406) 

where  the  "principal  factor"  T  is  a  function  of  the  form 


i 

with 


T 


(407) 


-  f  A  (X)p01 
i-1 


i Hx,p) 


(408) 


1 


w&mm 


the  q.  bain?  non-negative  decreasing  rational  numbers  and  p  being  a  positive  integer 
independent  of  j.  u  stands  for  tha  "secondary  factor" 

j.  . 

u(x,p)  -  l  *  (x)p  M  •  <409) 

v*0 

Here  M  is  a  positive  Integer  independent  of  j.  The  positive  integer  j'  depends  on 

j  and  increases  with  }.  (In  the  application  to  differential  equation  (101)  the  first 

term  of  (409)  is  obtained  if  j  »  -  - ,  as  we  shall  see.) 

n-« 

As  we  have  said  before,  the  maximal  value  of  j  for  which  a  solution  Y(x,p)  of 
(403)  can  be  constructed  depends  on  the  number  of  times  the  coefficients  Big  can  be 
differentiated.  It  can  be  only  determined  in  the  course  of  the  successive  construction  of 
the  terms  of  (409).  If  the  Big  can  be  differentiated  indefinitely,  then  j  and  j'  can 
have  arbitrarily  large  values. 

Remark i  The  theory  remains  valid  if  the  series  (40S)  are  not  convergent  but  only 
asymptotic  expansions,  nut  we  do  not  need  this  case  for  our  application. 
functional  Part:  bet  Y(x,p)  »  T«u  be  a  solution  of  (403).  Then  there  la  a  solution 
y(x,p)  of  the  exact  differential  equation  (401)  such  that 
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y(x,p)  -  T(u  +  p  M  [0])  (410) 

and  this  equation  can  be  formally  differentiated  at  least  n-1  times,  i.e.  it  can  be 
differentiated  treating  the  symbol  [0]  as  if  it  were  a  constant. 


63.  The  construction  of  the  asymptotic  solution  of  the  differential  equation: 
A)  The  principal  factor. 

The  first  step:  In  ^  substitute  an  expression  of  the  form 


(compare  (407),  (408),  (409)).  The  result  of  this  substitution  is  an  expression  of  the 


*  T1  \ 
l  f  (X)P  ♦  P  *[0] 
i-1 


4 

4 

*4 


i...  v-.  -  vV  ’$u. 


If 


*  -v>. 


a 

the  first  term  A^(x)P  1  in  (408)  in  such  a  way  that  the  coefficients  fj(x)  of  the  term 
of  highest  order  in  (411)  vanishes* 

Noaillon  gives  a  general  method  that  allows  us  to  determine  the  exponent  and  the 

function  A^x)  systematically  in  the  general  case.  But  our  application  being  of  such  a 
simple  type,  and  A^(x)  can  be  found,  in  our  case,  more  directly. 

Second  step:  Having  chosen  the  first  term  A^P  of  ♦  we  set 


P 


y 


where  is  the  denominator  of  the  rational  number  .  This  substitution  transforms  the 


expression 


into  an  expression 


w 


We  now  repeat  the  reasoning  of  the  first  step  with 


L2(y2> 


instead  of  *****>,  in  order  to 

y 


find  the  term  A  P  ,  considering,  of  course,  for  a  only  values  that  are  smaller  than 
2  2  a 

a.  we  continue  in  this  manner  until  we  arrive  at  the  last  term  A  P  P  for  which  a  is 
1  P  P 

still  positive.  That  there  is  such  a  term,  i.e.  that  we  always  attain  an  exponent 
a  .  <  0  after  a  finite  number  of  steps,  is  proved  in  Noail Ion's  paper. 

p+l 

B)  The  secondary  factor. 

First  step:  In  L(Y)  substitute  Y  -  T'u  where  T  is  the  function  calculated  in  A) 

and  u  as  yet  undetermined.  The  result  of  the  substitution  is  an  expression  of  the  form 

T  K(u,P),  where  K(u,P)  is  a  linear  differential  expression  in  u  whose  coefficients  are 
1  /M 

power  series  in  o  «  P  ,  M  being  the  common  denominator  of  all  the  exponents  in 

♦.  In  K(u,P)  collect  the  terms  of  highest  order  in  o.  Then  K(u,P)  can  be  written  in 


the  form 
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K(u,P)  -  S  G(u)  +  H(U/P) 


(412) 


Here  S  ia  the  highest  power  of  o  occurring  in  K(u,P)  and  G(u)  and  H(u,P)  are 
differential  expressions,  G(u)  being  independent  of  P.  It  can  be  proved  that  the 
highest  order  of  differentiation  occurring  in  G(u)  ia  greater  than  tero. 

Second  step:  Find  a  solution  of  the  differential  equation 

G(u)  -  0  .  (413) 

We  take  this  solution  as  the  first  term  Eq(x)  of  the  series  (409).  Since  we  want  Bg(x) 
to  be  bounded  in  the  whole  interval  of  x  in  which  we  consider  the  asymptotic  expansion, 
we  have  to  add  to  the  assumptions  1*  -  5»  in  section  62  the  condition  6*:  The  coefficient 
of  the  highest  derivative  in  (413)  does  not  vanish  in  any  part  of  the  interval  ot  <  x  <  6. 

Third  step:  In  order  to  find  the  function  E^(x),  W  >  0  of  (409),  Noaillon  proceeds 
as  follows:  He  determines  by  recursion  a  sequence  of  functions  w^(x,  )  w j(x,  ),...  from 
the  formula 

G(wv>  -  |  H(wv_i#P)  .  (414) 

It  is  easily  seen  that  the  w^tx.P)  are  of  the  form 

W  -  I  s"1  gv. (x)  .  (415) 

i-0 

The  functions  E^lx)  are  then  given  by 

P 

eu(x)  “  I  9VU(*>  .  (416) 

v*  i 

It  is  not  difficult  to  prove  that  the  function  (406),  if  determined  by  the 
construction  which  we  have  just  outlined  here,  does  in  fact  satisfy  the  relation  (403). 

The  construction  of  the  asymptotic  solution  Y(x,P)  is  by  no  means  uniquely 
determined,  and  it  can  be  proved  that  the  construction  yields  asymptotic  expansions  which 
are  asymptotic  approximations  to  a  fundamental  system  of  solutions  of  (401). 
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64.  Asymptotic  Solution  of  L(y)  =  N(y)  +  M(y)  -  0. 

A)  The  principal  factor. 

First  step:  We  now  apply  the  construction  described  in  section  63  to  the  particular 
type  of  differential  equation  under  consideration.  Substituting 

y  »  Y  »  T»u 


MXi£i(  where  L(y,p)  is  now  the  differential  expression  (101),  and  T  and  u  are 

expressions  of  the  form  described  by  (407),  (408),  (409),  we  see  that,  unless  i|i  «  0,  the 

term  of  highest  order  of  ^ -  is  (Aj(x)p  )  .  The  condition  that  the  highest  terms  of 

y  N(  y)  M(  y) 

- as  function  of  p,  and  of  — *—  cancel  out  is  therefore 

P  Y  Y 

a.n-1  a.m 

a"(x)p  h  -bQ(x)  a“(x)p  , 


o^n  -  1  -  a^ 


(417) 


and 


»n  V. 

1  “  "Vi  ' 


(418) 


provided  iji  /  0.  From  (417)  and  (418)  we  conclude 

1 


1  n-a 


(419) 


and 


A,  -  (-b0)1/n_n 


(420) 


Second  step:  Now  we  write 


o  ft 


y  =»  e 


(421) 


where  o  is  defined  by  p  3  a  and  (x)  is  one  of  the  functions  y(x)  defined  in 
theorem  1,  section  6.  Yj  is  the  function 


ft  MOdC 

0 

y2  -  e  •  u 


Here  u  is  again  the  series  (409)  and 

P 

<M*>  -  l  A  (x)  p 
i-2 


(422) 


(423) 
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He  find 


v(v)  „  oW-1^1(x)y- 

1 -  *  o  (x)  +  V  —  +  •  •  •  • 

y  y. 


(424) 


The  term  of  highest  order  in  the  second  term  of  the  right  member,  if  Y^  is  replaced  by 
its  value  according  to  (422)  and  (423)  has  the  value 

,  v-1  °2  v-1 

v  A  2  a  p  if 

V  V 

This  term  is  of  lower  order  in  p  than  the  term  a  <t  ,  since 

*(»-»)«, 

p  o  ■  o  -1 


and 


in  consequence  of  the  assumption 


( n-m)  ol  -  1  <  0  , 


^  <  °t  "  n-m  • 

The  omitted  tens  in  (424),  which  are  indicated  by  dots,  are  similarly  seen  to  be  of  still 

b(yj)  Inn  mm 

lower  order.  The  terms  of  highest  order  in  - ,  i.e.  —  a  and  b.o  «  cancel  out, 

y2  P  0 

since  i(x)  and  o  have  been  chosen  such  as  to  achieve  just  this.  The  next  terms  are 

(425) 


1  n-1  °2  n-1 

np  V  p  *  .  tor  -  y 


1  “‘V 


and 


„  „  -1  “2  -1  ,  M<y2> 

b A  o  p  <e  ,  for  — - 

u  *  y*i 


(426) 


provided  <jj  >  0.  (For,  if  “  0,  there  are  more  terms  of  the  same  order  as  (425)  and 
(426) .) 

Following  Noaillon's  construction  we  try  to  determine  and  A2  in  such  a  way  that 

these  two  terms  cancel  out.  But  setting  the  sum  of  (425)  and  (426)  equal  to  zero,  and 
inserting  for  *(x)  its  value  (-bg(x)  )1'/n“®  leads  to 

n-m  -  0  , 

which  was  excluded. 

Hence,  cij  >  0  is  impossible  and  therefore  •  <w(x). 
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B)  The  secondary  tern. 


He  follow  the  construction  of  Noaillon  In  order  to  find  the  differential  expression 
K(u,p)  of  (412).  Since 

(v)  v  v  .  vl  v-1  ,  . 

y  ■  Tlu  «ii  ♦  v«  v  u*  •»■•••  •) 

where  the  dots  Indicate  terns  of  lower  order,  we  find  by  an  easy  calculation 


L(Tu)  •  T  K(u,p)  “  T{om  ^(a^11  ^u  +  nf”  'u*  +  b^®  ^ u  +  nb^B”1u' )  +  H(u,p)} 


Hence,  the  function  Bg(x)  In  this  case  la  a  solution  of  the  differential  equation 
obtained  by  setting  the  factor  of  o  equal  to  zero.  This  differential  equation  can  be 

written 

(a^  ♦  b^Ju  *  (np  ♦  rab^u' 

or 

(a,  b0  -  b,)u  »  b0(n  -  m)u’  . 

Therefore 

fx  Vg)bo<C)  *  VS)  ir 

"  *o  b0<5)(n-n)  ^ 

B0(x)  -  n(x)  -  e  .  (427) 

Interested  In  the  other  terras  I^x),  v  >  0  of  (409). 

1*  -  6*  of  Noalllon's  theory  are  satisfied  In  our  application  for  the  whole 
<  8.  Condition  5»,  in  particular,  is  equivalent  to  condition  6»  of  the 
chapter  I.  Hence,  we  conclude  that  there  are  n-m  solutions  of  (101)  of 

o  F  *  <«« 

Yv(x,p)  -  e  a  v  (n(x)l,  *(v  -  1,2,...,n-m)  .  (42S) 

Note  that  the  function  n(x)  Is  the  same  for  all  Yy(x,  )  and  that  It  does  not  vanish  In 
a  <  x  <  6. 


He  are  not 
Conditions 
interval  a  <  x 
Main  Theorem  of 
the  form 
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But  w«  can  find  more  asymptotic  solution*  of  th«  dlffarsntlal  aquation  (101)  by 
dropping  th*  assumption  that  th*  i(i(x,p)  of  (408)  is  not  ssro.  In  fact.  If  th#  principal 
factor  T  of  (40?)  is  equal  to  1,  the  method  used  for  the  construction  of  th*  secondary 
factor  In  section  63,  B)  leads  to  asymptotic  solutions  given  by  series  of  the  fora  (409). 
The  first  term  of  each  of  these  series  is  a  solution  of  the  differential  equation 
M(y)  -  0.  Taking  a  fundamental  system  u^(x)  of  m  Independent  solutions  of  this 
differential  equation  as  first  terms  of  m  asymptotic  solutions  of  (101)  we  can  add  to  the 
n-m  solutions  of  (101)  given  by  (428)  a  more  solutions  of  the  form 

V(x,p)  -  lu  (x)}  ,  (v-  1,2, .. . ,m)  .  (429) 

n-a+v 

linerly  Independent,  for  sufficiently  large  p,  can  be  given  by  calculating  the  asymptotic 
value  of  the  Wronskian  of  these  n  functions.  It  does  not  offer  any  difficulties. 

This  finishes  the  proof  of  theorem  1  of  section  6. 
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